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Abstract 

The  Spectrally  Modulated,  Spectrally  Encoded  (SMSE)  framework  provides  an  ef¬ 
fective  means  for  implementing  Orthogonal  Erequency  Division  Multiplexing  (OEDM) 
signals  -  a  bedrock  technology  for  future  fourth  generation  (4G)  communication  systems 
based  on  Cognitive  Radio  (CR)  and  Software  Defined  Radio  (SDR)  techniques.  As  4G 
SMSE  communications  emerge  they  must  coexist  with  other  systems  while  competing  for 
available  communication  resources.  Given  a  lack  of  inter-system  orthogonality  and  limited 
available  bandwidth,  these  signals  must  be  designed  to  spectrally  coexist  while  inducing 
“manageable”  levels  of  mutual  interference.  The  research  goal  was  to  demonstrate  a  struc¬ 
tured  means  for  SMSE  waveform  design  using  two  techniques  commonly  employed  in 
operations  research:  Genetic  Algorithm  (GA)  and  Response  Surface  Methodology  (RSM). 

The  design  process  is  demonstrated  herein  for  a  coexistent  scenario  containing  SMSE 
and  Direct  Sequence  Spread  Spectrum  (DSSS)  signals.  Coexistent  SMSE-DSSS  designs 
are  addressed  under  both  perfect  and  imperfect  DSSS  code  tracking  conditions  using  a  non¬ 
coherent  delay-lock  loop  (DEE).  Under  both  conditions,  the  number  of  SMSE  subcarriers 
N f  and  subcarrier  spacing  A/  are  the  optimization  variables  of  interest.  Eor  perfect  DEE 
code  tracking  conditions,  the  GA  and  RSM  optimization  processes  are  considered  inde¬ 
pendently  with  the  objective  function  being  end-to-end  DSSS  bit  error  rate  P^-  A  hybrid 
GA-RSM  optimization  process  is  used  under  more  realistic  imperfect  DEE  code  tracking 
conditions.  In  this  case,  optimization  is  accomplished  through  a  correlation  degradation 
metric  with  the  GA  process  being  first  applied  to  generate  a  “coarse”  solution  followed  by 
RSM  processing  which  provides  the  final  optimized  solution. 

Eor  all  perfect  and  imperfect  DEE  code  tracking  scenarios  considered,  the  optimized 
DSSS  Pb  minimization  results  yielded  SMSE  waveform  designs  and  Pb  performance  that 
was  consistent  with  scenarios  having  no  coexistent  SMSE  signal  present  (best-case  co¬ 
existent  performance).  Eor  the  optimized  DSSS  Pb  maximization  solutions,  worst-case 
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SMSE-DSSS  coexistence  was  achieved  for  SMSE  waveform  designs  that  were  spectrally 
“matched”  to  the  DSSS  signal,  i.e.,  greatest  Pf,  degradation  was  experienced  when  the  re¬ 
sultant  SMSE  subcarrier  spacing  A/  was  an  integer  multiple  of  the  spectral  line  spacing 
A/chip  of  the  DSSS  spreading  code. 

This  work  has  successfully  expanded  the  practical  utility  of  a  previously  developed 
tool,  the  original  SMSE  framework,  by  demonstrating  a  more  efficient,  structured  means  for 
coexistent  waveform  design  that  replaces  previous  trial  and  error  methods.  The  research 
objective  has  been  achieved  in  the  sense  that  4G  communication  design  engineers  now 
have  one  additional  tool  at  their  disposal  and  its  significance  has  been  acknowledged  - 
the  technical  community  is  one  step  closer  to  actually  hitting  the  bedrock  of  OEDM-based 
signaling  using  the  SMSE  framework.  It  is  also  important  to  note  that  the  particular  DEE 
implementation  used  here,  and  the  metrics  used  to  characterize  various  tracking  conditions 
(perfect  and  imperfect),  are  sufficiently  general  such  that  the  optimization  demonstrations 
herein  are  broadly  applicable  to  other  non-communication  applications  employing  DEE 
tracking,  e.g.,  precision  navigation,  timing,  geolocation,  etc. 
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Application  of  Optimization  Techniques  to 
Spectraeey  Modueated,  Spectraeey  Encoded 
Waveeorm  Design 


I.  Introduction 

A  quick  glance  at  the  US  frequency  allocation  chart  [33]  shows  that  the  radio  spec¬ 
trum  between  9  KHz  and  300  GHz  is  fully  allocated.  However,  this  provides  no  indieation 
as  to  how  efficiently  the  allocated  speetrum  is  being  utilized.  As  the  need  for  inereasing 
data  rates,  the  number  of  users,  and  overall  network  eapaeity  improvement  has  grown,  the 
issue  of  speetrum  effieieney  versus  speetrum  alloeation  has  gained  mueh  attention.  This 
is  highlighted  by  the  following  relevant  sampling  of  reeent  literature  addressing  speetrum 
effieieney: 

1.  [39]  indicates  there  is  an  “underutilization  of  the  radio  speetrum  as  revealed  by 
extensive  measurements  of  aetual  speetrum  usage”  in  [12]. 

2.  [38]  eites  measurements  in  [12]  as  well  and  conclude  that  “at  any  given  time  mueh 
of  the  prized  speetrum  lies  idle  ...  speetrum  shortage  results  from  the  spectrum  man¬ 
agement  poliey  rather  than  the  physical  scarcity  of  usable  frequeneies.” 

3.  [9]  eites  measurements  in  [12]  as  well  and  indieate  that  “...  at  any  time  roughly  10% 
of  the  unlieensed  frequeney  spectrum  is  actively  in  use  (leaving  90%  unused)” 

4.  [8]  provides  an  overview  of  the  Defense  Advaneed  Researeh  Projeets  Ageney  (DARPA) 
Next  Generation  (XG)  Communieation  Program  and  indieates  that  94%  of  the  spee¬ 
trum  was  unused  worldwide  at  the  time  of  a  2002  study. 

5.  [1]  indicates  that  “temporal  and  geographical  variations  in  the  utilization  of  assigned 
speetrum  ranges  from  15%  to  85%  with  a  high  varianee  in  time”  and  eite  [25]  in 
support  of  their  eonelusion. 
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While  there  may  be  some  uneertainty  of  the  eurrent  value  of  speetrum  effieieney,  as  well 
as  its  terminology  (pereentage  used,  unused,  underused,  underutilized,  ete.),  the  overall 
eonsensus  is  that  speetrum  effieieney  is  mueh  poorer  than  what  ean  be  aehieved,  and  the 
teehnieal  eommunity  is  unified  in  its  quest  to  improve  it. 

To  some  degree,  the  teehnieal  eommunity  has  unified  under  the  eoneept  of  Dynamie 
Speetrum  Aecess  (DSA)  which  is  highlighted  through  recent  events  such  as  the  IEEE’s 
Dynamic  Spectrum  Access  Networks  (DySPAN)  symposiums  [10].  DySPAN  has  grown  to 
be  “the  preeminent  event  to  gather  international  economists,  engineers,  network  architects, 
researchers  and  academic  scholars  together  to  share  cutting  edge  research  on  and  demon¬ 
strations  of  emerging  wireless  technology.”  The  original  DySPAN  symposium  was  held  in 
2005  and  activities  since  then  have  continued  to  have  significant  international  influence  on 
policy,  technology,  research  and  development  of  next  generation  wireless  systems.  Prom  a 
technology  focused  perspective,  the  ability  to  effectively  employ  DSA  techniques  is  highly 
linked  to  Software  Defined  Radio  (SDR)  and  Cognitive  Radio  (CR)  concepts. 

1.1  Research  Motivation 

The  fundamental  DSA  concept  provides  one  means  for  fourth  generation  (4G)  and 
subsequent  communication  systems  to  improve  spectrum  efficiency  while  minimizing  ad¬ 
verse  coexistence  effects.  This  is  done  by  monitoring  and  adapting  to  changing  channel 
conditions,  traffic  requirements,  coexisting  signals,  and  spectrum  availability  by  generat¬ 
ing  waveforms  that  dynamically  respond  to  these  conditions.  The  coding,  modulation, 
and  multiple  access  techniques  of  emerging  systems  will  require  this  type  of  adaptivity,  as 
enabled  at  the  physical  layer  through  SDR  techniques.  The  SDR  advantages  are  obvious 
given  that  the  radio  can  be  easily  upgraded  with  changes  in  standards,  frequency  allocation, 
security  practices  and  real-time  environmental  changes.  To  some  degree,  these  advantages 
have  been  successfully  exploited  in  some  systems  [13, 19,23,24].  However,  these  systems 
are  somewhat  restricted  in  the  decision  stage  of  adaptation,  i.e.,  the  “brain”  that  controls 
the  adaptation  based  on  observed  environmental  factors  is  somewhat  limited. 
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To  address  this  limitation,  CR  techniques  have  been  proposed  and  investigated  to  in¬ 
crease  a  given  radio’s  autonomy  by  increasing  its  ability  to  observe,  “think”,  and  decide 
the  next  best  course  of  action.  This  represents  one  form  of  “brain  empowered  communica¬ 
tions”  referred  to  by  Haykin  [17].  In  this  context,  CR  principles  are  considered  a  means  for 
enhancing  and  advancing  SDR  functionality  and  capability.  By  way  of  maintaining  consis¬ 
tency  with  previous  work  in  [26, 30]  that  forms  the  basis  for  this  research,  the  synergistic 
union  between  CR  and  SDR  will  be  referred  to  here  as  CR-based  SDR.  In  this  context,  the 
SDR  is  a  software  controlled  waveform  generator  that  is  controlled  by  CR  principles  de¬ 
signed  to  improve  environmental  assessment  and  spectrum  usage  efficiency.  In  the  context 
of  achieving  smarter  communications,  Haykin  identifies  Orthogonal  Frequency  Division 
Multiplexing  (OFDM)  as  a  bedrock  technology  for  CR-based  SDR  implementation  [17]; 
as  with  previous  work  in  [26,30],  this  continues  to  provide  the  motivation  for  the  research 
presented  here. 

While  the  concept  of  OFDM  is  relatively  simple,  the  ability  to  generate  OFDM  sig¬ 
nals  has  only  recently  become  practical  from  a  hardware  perspective.  The  speed  of  modern 
digital  processors  now  allows  waveform  characteristics  to  be  defined  in  the  frequency  do¬ 
main,  with  conversion  to  time  domain  waveforms  via  an  Inverse  Fast  Fourier  Transform 
(IFFT)  occurring  at  rates  that  enable  high-speed  communications.  As  such,  OFDM  has 
proven  to  be  very  flexible  and  has  generated  significant  interest  throughout  the  research 
community  [11,20,34-37,40].  The  rapidly  expanding  pool  of  emerging  OFDM-based 
techniques,  as  illustrated  in  the  cloud  region  of  Fig.  1.1,  drove  the  need  to  develop  a  unified 
framework  to  encapsulate  OFDM  variants.  Elements  of  the  resultant  framework  are  cap¬ 
tured  in  the  analytic  expression  in  Fig.  1.1  which  effectively  embodies  a  class  of  what  are 
now  called  Spectrally  Modulated,  Spectrally  Encoded  (SMSE)  waveforms  [26-30]. 

Researchers  continue  to  investigate  expanded  roles  for  applying  the  original  SMSE 
framework.  Some  of  the  more  recent  efforts  are  focusing  on  overlay,  underlay  and  hybrid 
overlay /underlay  waveform  implementations  [5,6].  In  the  context  of  an  overlay  wave¬ 
form  architecture  where  interference  to  primary  users  is  mitigated  by  avoiding  their  spec¬ 
tral  regions,  these  latest  SMSE  developments  suggest  that  the  original  SMSE  framework 
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Unifying  SMSE  Framework 


Figure  1.1:  Unifying  SMSE  framework  for  OFDM-based  signaling  using 
a  CR-based  SDR  architeeture.  Waveform  adaptivity  is  provided  through 
variation  in  design  variables  denoted  by  [  d  c  w  o  a  u  ]  [26]. 


is  best  eharaeterized  as  applying  hard  decision  eriteria  to  speetral  usage.  By  accounting 
for  possible  implementation  of  both  overlay  and  underlay  techniques,  these  recent  activi¬ 
ties  have  demonstrated  preliminary  success  with  what  is  being  called  soft  decision  SMSE 
(SD-SMSE).  While  the  research  focus  in  this  dissertation  is  on  the  original  SMSE  frame¬ 
work  and  its  optimization  for  coexistence  with  other  systems,  emergence  of  SD-SMSE 
techniques  certainly  warrants  future  investigation  given  that  the  methods  considered  here 
should  to  be  directly  applicable. 


1.2  Research  Assumptions 

Definitions  of  mathematical  symbols  and  terminology  are  provided  throughout  the 
document  as  they  are  introduced.  Results  of  all  work  presented  should  be  interpreted  within 
limits  and  constraints  imposed  by  the  assumptions  that  have  been  made.  The  following 
summarizes  the  assumptions  that  made  for  the  research. 

1.2.1  Coexistent  SMSE-DSSS  Scenario. 

•  For  all  scenarios  considered,  the  coexistent  SMSE  and  DSSS  signals  are  spectrally 
coexistent  (same  center  frequency)  and  are  operating  over  an  Additive  White  Gaus¬ 
sian  Noise  (AWGN)  channel.  However,  given  the  focus  of  the  waveform  design 
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demonstrations  are  on  process,  the  procedures  used  are  valid  for  other  channel  mod¬ 
els. 

•  For  coexistent  and  interference  scenarios,  the  relative  power  ratios  between  the  signal 
of  interest,  the  interfering  signal,  and  the  channel  noise  are  the  dominant  factors. 
Therefore,  the  research  assumes  these  relative  power  levels  are  set  such  that  SMSE 
parametric  changes  within  can  impact  DSSS  receiver  performance. 

•  Except  for  spreading  code  tracking,  the  DSSS  receiver  was  perfectly  synchronized 
to  the  transmitted  DSSS  signal  in  terms  of  carrier  tracking  (frequency  and  phase) 
and  communication  symbol  tracking.  In  addition,  final  results  presented  in  Sec¬ 
tion  4.1  are  based  on  perfect  spreading  code  tracking  while  the  results  presented 
in  Section  4.2  are  based  on  imperfect  spreading  code  tracking. 

1.2.2  Optimization  Techniques. 

•  Eess-than-global  optimal  solutions  were  deemed  acceptable  for  demonstration  pur¬ 
poses.  The  optimal  responses  to  SMSE  variable  changes  were  at  least  locally  optimal 
and  accepted  given  the  research  goal  was  to  find  “good”  parameter  values  without  re¬ 
quiring  exhaustive  testing. 

•  All  GA  experiments  were  designed  to  ensure  that  all  possible  combinations  of  SMSE 
parameters  {N f  ,A/)  were  allowed,  i.e.,  all  combinations  were  in  the  optimization 
feasibility  region. 

•  All  RSM  experiments  used  a  second-order  /3  model  for  optimization.  The  experi¬ 
ments  consisted  of  a  two-factor,  three-level,  full-factorial  design  with  four  additional 
center  runs.  The  SMSE  parameters  were  optimized  using  the  steepest  ascent/descent 
process  until  the  response  surface  fit  the  second-order  model,  as  determined  by  an 
ANOVA. 
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1.3  Research  Sponsorship 

This  research  was  sponsored  in  part  by  the  Sensors  Directorate  of  the  Air  Force  Re¬ 
search  Laboratory  (AFRL/RY),  Wright-Patterson  AFB,  Ohio.  The  work  performed  and 
results  obtained  directly  support  their  vision  of  providing  sensor  and  countermeasure  tech¬ 
nology  to  enable  complete  freedom  of  air  and  space  operations  for  the  military  warfighter 
and  for  civilian  agencies  supporting  homeland  security.  The  work  most  closely  aligns  with 
the  discovery  and  development  aspects  of  AFRL/RY’s  mission  which  is  targeted  toward 
producing  affordable  sensor  and  countermeasure  technologies. 
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II.  Background 


2.1  Digital  Communications 

Electronic  communications  can  occur  using  either  analog  or  digital  waveforms.  Ana¬ 
log  communication  applications  include  television  (TV),  AM  radio,  FM  radio,  and  early 
cellular  telephones.  Digital  communication  techniques  differ  from  analog  techniques  in 
that  information  is  sent  using  a  set  of  predefined  quantized  values  that  are  commonly  rep¬ 
resented  by  either  a  binary  value  of  0  or  1 .  Digital  communication  applications  include 
modern  cellular  telephone  systems,  high  definition  television  (HDTV)  and  a  multitude 
of  wireless  computer  network  implementations  [31].  There  are  many  reasons  why  mod¬ 
em  communication  system  designers  prefer  digital  techniques.  First,  the  digital  receiver 
must  only  decide  between  a  finite  number  of  transmit  conditions  (communication  sym¬ 
bols).  Thus,  digital  techniques  perform  very  well  in  noisy  environments  since  the  original 
signal  does  not  have  to  be  reconstructed.  Second,  digital  repeater  systems  may  be  imple¬ 
mented  as  regenerative  receive-transmit  nodes.  They  are  able  to  receive,  estimate  symbols, 
perform  error  correction  and  duplicate  information  before  passing  it  on,  allowing  reliable 
long  distance  communications.  Also,  digital  hardware  is  often  less  expensive  to  manufac¬ 
ture  than  analog  equipment  [31]. 

A  basic  digital  communications  system  model  is  shown  in  Fig.  2.1.  The  digital  input 
sequence  {dk},  dk  G  {0, 1}  ,  is  mapped  to  a  waveform  s{t)  for  transmission.  This  process 
of  mapping  information  to  a  transmitted  waveform  is  called  modulation.  The  transmitted 
signal  propagates  through  the  channel,  or  transmission  medium,  and  is  corrupted  by  noise. 
The  channel  may  be  a  transmission  line,  such  a  telephone  line  or  coaxial  cable,  or  the 
open  airwaves  as  assumed  for  this  research.  The  term  noise  describes  the  combined  effect 
of  interference  from  many  sources,  including  the  combined  effect  of  additional  signals  in 
the  channel,  atmospheric  effects  and  thermal  effects  within  the  receiver  itself.  Additive 
white  Gaussian  noise  (AWGN)  is  a  common  channel  model  used  for  designing  and  ana¬ 
lyzing  electronic  communications  systems.  As  used  for  this  research,  the  AWGN  channel 
assumption  provides  a  reasonable  starting  point  for  system  design. 
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Figure  2.1:  Basie  digital  eommunieation  system  model  showing  trans¬ 
mitter  modulation  and  reeeiver  demodulation  funetions.  An  AWGN  propa¬ 
gation  ehannel  is  illustrated. 

M-ary  Phase  Shift  Keying  (MPSK)  is  one  type  of  digital  data  modulation  whereby 
the  input  data  bits  are  mapped  to  a  series  of  transmitted  communication  symbols  given  by 


s{t)  =  A  cos  [271  fct  +  (p  (f)]  , 


(2.1) 


where  0  <  f  <  Tgym,  Tgym  is  the  symbol  duration,  and  phase  value  (j){t)  is  determined  by 
the  input  bit  values.  The  special  case  where  each  symbol  represents  M  =  2  bits,  is  known 
as  Binary  Phase  Shift  Keying  (BPSK)  which  is  the  basis  for  the  analysis  in  this  dissertation. 
In  this  case,  0  =  0°  {dk  =  1)  or  <pi  =  180°  {dk  =  0)  and  the  expression  in  (2.1)  can  be 
rewritten  as 


s{t)  =  (  —  1)'^'“  A  cos  {271  fct)  . 


(2.2) 


The  plot  in  Fig.  2.2  shows  a  representative  BPSK  waveform  spanning  two  symbol  periods 
having  different  modulation  values  (either  dk  =  0  or  dk  =  1).  This  bit  difference  causes 
the  180°  phase  shift  occurring  at  Tgym  in  the  plot. 

To  use  the  available  communication  resources  efficiently,  digital  waveforms  are  of¬ 
ten  multiplexed  within  or  across  specific  signaling  domains  (time,  frequency,  space,  po¬ 
larization  and/or  code).  Several  multiplexing  schemes  exist  that  allow  multiple  users  to 
effectively  share  communication  resources,  including  Time  Division  (TDM),  Frequency 
Division  (FDM),  Space  Division  (SDM),  Polarization  Division  (PDM)  and  Code  Division 
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Figure  2.2:  Representative  Binary  Phase  Shift  Keyed  (BPSK)  waveform 
spanning  to  symbol  periods  having  different  data  modulation  values.  The 
instantaneous  180°  phase  change  at  Tsy^  is  due  to  this  difference  [31]. 

(CDM)  [31].  Of  particular  interest  to  this  research  are  FDM  and  CDM,  each  of  which  is 
introduced  and  discussed  in  greater  detail  where  appropriate. 

2.1.1  Spectrally  Modulated,  Spectrally  Encoded  (SMSE)  Framework.  One  method 
to  allow  more  data  through  a  given  channel  involves  dividing  the  allocated  frequency  band 
into  several  narrower  subbands,  allowing  multiple  individual  signals  to  coexist.  This  tech¬ 
nique  is  known  as  Frequency  Division  Multiplexing  (FDM)  [31].  A  representative  FDM 
spectrum  is  illustrated  in  Fig.  2.3  for  three  subbands  separated  by  A/.  The  receiver  for 
a  particular  signal  in  this  scheme  is  able  to  extract  the  desired  signal  using  an  appropri¬ 
ate  band  pass  filter  centered  on  the  subband  of  interest  and  having  a  bandwidth  that  is  a 
fraction  of  A/,  say  Wbp  =  A//2.  In  this  case,  a  “manageable”  amount  of  interfering 
signal  power  from  the  other  two  subbands  passes  through  the  filter  and  impacts  demodula¬ 
tion  [31].  Given  that  the  signals  are  separated  in  frequency,  the  important  design  parameters 
for  FDM  include  frequency  separation  A/  and  Wbp- 

One  particularly  efficient  way  to  space  FDM  subcarrier  frequencies  is  to  assign  A/ 
such  that  the  individual  subband  center  frequencies  are  spectrally  coincident  with  the  first 
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null  of  adjacent  signal  [34].  For  this  to  occur,  the  symbol  duration  and  subearrier  separation 
are  related  such  that  Tsym  =  l/^/>  i-e->  the  individual  subearrier  frequeneies  are  harmon¬ 
ies  of  the  fundamental  frequeney  A/.  This  partieular  subcarrier  spaeing  is  illustrated  in 
Fig.  2.4.  Ideally,  this  proeess  results  in  individual  signals  being  mathematieally  orthogonal 
with  no  mutual  interferenee.  This  type  of  FDM  seheme  is  known  as  Orthogonal  Frequeney 
Division  Multiplexing  (OFDM)  [34] . 

While  the  concept  of  OFDM  is  relatively  simple,  the  ability  to  generate  OFDM  sig¬ 
nals  has  only  reeently  beeome  praetieal  from  a  hardware  perspeetive.  The  speed  of  modern 
digital  proeessors  now  allows  waveform  eharaeteristies  to  be  defined  in  the  frequeney  do¬ 
main,  with  conversion  to  time  domain  waveforms  via  an  Inverse  Fast  Fourier  Transform 
(IFFT)  occurring  at  speeds  that  enable  high-speed  communications.  As  such,  OFDM  has 
proven  to  be  very  flexible  and  has  generated  significant  interest  throughout  the  researeh 
community  [11,20,34-37,40].  The  rapidly  expanding  pool  of  emerging  OFDM-based 
teehniques  drove  the  need  to  develop  a  unified  framework  to  eneapsulate  OFDM  variants. 
The  resultant  framework  effeetively  embodies  a  elass  of  what  are  now  ealled  Speetrally 
Modulated,  Speetrally  Eneoded  (SMSE)  waveforms  [27-30]. 

Researchers  continue  to  investigate  expanded  roles  for  applying  the  original  SMSE 
framework,  with  some  of  the  more  reeent  efforts  foeusing  on  hybrid  overlay /underlay  wave¬ 
form  implementations  [5,  6].  In  the  eontext  of  an  overlay  waveform  arehiteeture  where 
interferenee  to  primary  users  is  mitigated  by  avoiding  their  speetral  regions,  these  latest 
SMSE  developments  suggest  that  the  original  SMSE  framework  is  best  eharaeterized  as 
applying  hard  decision  eriteria  to  spectral  usage.  By  aeeounting  for  possible  implemen¬ 
tation  of  both  overlay  and  underlay  techniques,  these  reeent  activities  have  demonstrated 
preliminary  sueeess  with  what  is  being  ealled  soft  decision  SMSE  (SD-SMSE).  While  the 
researeh  foeus  in  this  dissertation  is  on  the  original  SMSE  framework  and  its  optimization 
for  eoexistenee  with  other  systems,  emergenee  of  SD-SMSE  teehniques  eertainly  warrants 
future  investigation  given  that  the  methods  eonsidered  here  should  be  direetly  applieable. 
Eor  eompleteness,  the  following  SMSE  development  is  provided  and  is  based  on  the  orig- 
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Figure  2.3:  Representative  Frequency  Division  Multiplexed  (FDM)  spec¬ 
trum  for  three  subcarriers  separated  in  frequency  by  A/  [34]. 


Frequency 

Figure  2.4:  Representative  SMSE  signal  spectrum  for  three  subcarriers 
separated  in  frequency  by  A/  =  1/Tsym^  an  integer  multiple  of  the  symbol 
interval  [34]. 
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inal  SMSE  framework  presented  in  [30].  The  reader  is  referred  to  the  original  SMSE 
framework  documentation  if  additional  details  are  required. 

Accounting  for  all  SMSE  waveform  design  variables,  the  original  SMSE  framework 
provides  a  unified  approach  for  generating  and  characterizing  a  host  of  OEDM-based  sig¬ 
nals.  Using  0  to  denote  Hadamard  product  (element-by-element  multiplication),  the  spec¬ 
tral  representation  of  the  SMSE  symbol  is  given  by  [28, 30] 


Sfc  =  c  0  dfc  0  w  0  Ofc  ,  (2.3) 

where  the  design  variables  are  defined  as  follows:  Coding  c  =  [ci,  C2, . . . ,  catj.],  Q  G  C; 
Data  Modulation  d  =  [di,  d2,  ■  ■  ■ ,  d^^],  di  G  C;  Windowing  w  =  [tui,  W2,  ■  ■  ■ ,  WN^],Wi  G 
C;  Orthogonality  o  =  [oi,  02,  •  •  • ,  onj],  o*  G  C,  |oi|  =  1  V  i.  Each  of  these  terms  are 
introduced  to  functionally  incorporate  various  waveform  design  characteristics  commonly 
employed  in  communications. 

All  that  remains  to  completely  specify  the  SMSE  waveform  is  the  frequency  compo¬ 
nent  selection  and  symbol  duration  of  the  resulting  waveform.  The  frequency  component 
defines  the  number  of  carrier  components  that  are  spectrally  modulated  and  encoded.  As¬ 
suming  an  -point  inverse  fast  Eourier  transform  (lEET),  there  are  initially  Nf  possible 
carrier  components  available.  Use  of  components  from  this  pool  of  frequencies  is  con¬ 
trolled  through  frequency  assignment  and  use  variables.  Eor  example,  a  system  may  elect 
to  assign  a  subset  of  Nf  carriers  to  a  given  user.  This  assignment  of  frequencies  is  ac¬ 
counted  for  through  variable  a  =  [oi,  02, . . . ,  a*  G  {0, 1},  where  zeros  indicate  unas¬ 
signed  carriers.  Erom  this  assigned  pool  of  carriers,  some  may  be  unused  due  to  excessive 
interference,  system  design,  etc.  The  remaining  used  carriers  are  accounted  for  through 
variable  u  =  [ui,  M2, . . . ,  mat/],  w*  G  {0, 1},  where  zeros  indicate  unused  carriers  and  there 
are  total  of  P  used  carriers  (P  <  N f).  Thus,  u  is  a  subset  of  a  (u  C  a)  and  only  assigned 
carrier  components  exist  SMSE  symbols.  The  frequency  assignment  and  use  variables  are 
incorporated  into  (2.3)  as  follows 
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Sfc  =  a0u0c0dfc0w©Ofc, 


(2.4) 


where  the  carrier  component  of  is  given  by 

Sk[m]  =  .  (2.5) 

There  are  m  =  0,1, ...  ,Nf  -  1  frequency  components  with  c^,  Oc^,  dyn,k,  Wm, 

Ou,^  and  9o^  ^  being  the  corresponding  magnitudes  and  phases  of  the  design  variables.  As 
indicated  by  the  subscripted  k  and  m  indices  in  the  righthand  side  of  (2.5),  the  coding 
and  windowing  terms  only  vary  with  frequency  index  m,  whereas  the  data  modulation  and 
orthogonality  terms  vary  with  symbol  index  k  as  well. 

The  block  diagram  in  Fig.  2.5  illustrates  the  functional  processes  that  are  commonly 
used  to  generate  OFDM  signals.  Input  data  bits  dk  €  {0, 1}  are  mapped  to  discrete  BPSK 
coefficients  according  to  6^  =  (—1)'^'=.  The  BPSK  coefficients  are  first  grouped  using  a 
serial-to-parallel  (S/P)  conversion  process  and  then  weighted  by  vector  Wk  prior  to  the  IFFT 
operation.  Following  the  IFFT  operation,  the  grouped  bits  are  converted  back  to  a  stream 
stream  via  parallel-to-serial  (P/S)  conversion.  The  resultant  complex  discrete  samples  are 
then  converted  to  in-phase  and  quadrature  baseband  signals  by  digital-to-analog  conversion 
(DAC).  Finally,  the  complex  components  are  combined  in  a  quadrature  modulator  and  up- 
converted  to  the  desired  carrier  frequency  fc  for  transmission.  For  the  coexistence  analysis 
in  this  work,  the  key  SMSE  waveform  design  variables  include  the  number  of  IFFT  points 
Nf,  the  subcarrier  frequency  spacing  A/  and  the  inter-subcarrier  complex  weighting  Wk. 

2.1.2  Direct  Sequence  Spread  Spectrum  (DSSS)  System.  Spread  spectrum  com¬ 
munications  are  a  class  of  signals  that  employ  Code  Division  Multiple  Access  (CDMA). 
As  a  mechanism  for  enabling  multiple  access,  CDMA  provides  the  ability  to  increase  the 
number  of  users  within  a  finite  allocation  of  spectrum.  In  terms  of  multiple  access  perfor¬ 
mance,  CDMA  has  an  inherent  advantage  over  Time  Division  Multiple  Access  (TDMA) 
given  there  is  no  need  for  precise  timing  between  users  in  the  network.  Two  other  desirable 
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Figure  2.5:  Block  diagram  illustrating  functional  processes  commonly 
used  for  OFDM  symbol  generation  and  transmission  [34] . 

characteristics  of  CDMA  when  implemented  with  spread  spectrum  techniques,  include: 

1)  users  enjoy  some  degree  of  privacy  from  unintended  eavesdroppers  given  that  the  code  is 
only  shared  with  intended  receivers,  2)  spreading  of  information  across  a  wide  band  of  fre¬ 
quencies  provides  some  inherent  protection  against  fading  channel  effects,  and  3)  the  initial 
despreading  operation  in  the  DSSS  receiver  provides  some  level  of  protection  against  coex¬ 
isting  interferers  or  jammers  [25,31].  Spectral  spreading  in  a  DSSS  is  generally  achieved 
using  a  psuedorandom  noise  (PN)  code.  The  PN  code  is  actually  deterministic  but  pos¬ 
sesses  some  properties  that  are  similar  to  random  noise  when  observed  without  knowledge 
of  the  code  [31]. 

A  typical  DSSS  waveform  generation  process  is  illustrated  in  Fig.  2.6.  For  BPSK 
data  modulation,  the  generation  process  first  begins  by  mapping  the  sequence  of  input  data 
bits  {dk}  to  waveform  values  to  created  the  data  modulated  signal  d{t)  given  by 

-l)‘^>=pit)  ,  (2.6) 

/  G  (0,  Tgym] 

otherwise 

The  resultant  data  modulated  waveform  d{t)  modulates  the  carrier  signal  A  cos  {271  fct), 
where  fc  is  the  carrier  frequency.  The  carrier  modulated  signal  is  then  spectrally  spread  by 


where 


d{t)  = 


p{t)  = 


0, 
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Figure  2.6:  Typical  DSSS  signal  generation  process.  The  carrier  carrier 
fc  is  first  modulated  by  data  dependent  waveform  d{t)  and  then  spectrally 
spread  by  coded  waveform  c{t)  [25]. 

the  PN  coded  waveform  c{t).  Although  not  a  requirement,  the  spreading  waveform  c{t)  is 
generally  created  using  the  same  BPSK  technique  as  used  for  d{t).  Thus,  given  a  binary 
PN  spreading  sequence  of  elements,  {c}  =  {ci,C2,  one  period  of  spreading 

waveform  c{t)  contains  Nc  chip  intervals  of  duration  Tc  and  is  given  by 


c{t)  =  {-iT'^v{t) , 


(2.7) 


where 


1,  t  E  (0,  Tc] 
0,  otherwise 


v{t) 


The  plots  in  Fig.  2.7  graphically  illustrate  the  DSSS  waveform  construction  process. 
The  plots  in  Fig.  2.7a  and  Fig.  2.7b  represent  the  baseband  data  modulated  and  PN  coded 
signals  d{t)  and  c{t),  respectively.  These  signals  were  generated  using  the  BPSK  bit-to- 
waveform  mapping  processes  detailed  in  (2.6)  and  (2.7).  For  visual  clarity,  only  four  chip 
intervals  per  symbol  duration  were  used  {Tsym  =  1/Rsym  =  4  x  Tc  =  4/i?c).  The  signal 
in  Fig.  2.7c  is  the  result  of  multiplying  the  data  modulated  signal  d{t)  with  the  RF  carrier 
A  cos  {271  fct).  For  visual  clarity,  only  eight  carrier  frequency  cycles  per  symbol  duration 
were  used  (Tsym  =  8/  fc).  The  carrier  phase  transitions  in  modulated  carrier  are  read¬ 
ily  apparent  and  correspond  directly  with  phase  transitions  in  d{t).  Finally,  the  signal  in 
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Figure  2.1  \  Waveform  construction  process  for  BPSK  DSSS  system: 

(a)  baseband  data  modulated  waveform  d{t),  (b)  baseband  spreading  modu¬ 
lation  waveform  c(f),  (c)  carrier  modulated  waveform  d{t)  cos  {271  fj:),  and 
(d)  final  DSSS  waveform  c{t)d{t)  cos  {27r f^t)  [25]. 

Fig.  2.7c  is  modulated  by  the  PN  coded  signal  c{t)  to  create  the  transmitted  signal  shown 
in  Fig.  2.7d. 

The  spectral  spreading  resulting  from  this  process  is  illustrated  in  Fig.  2.8  which 
shows  power  spectral  density  (PSD)  plots  for  the  unspread  data  modulated  carrier  (dashed 
line  with  peak  value  Pbpsk)  and  the  final  spread  waveform  (solid  line  with  peak  value 
Pdsss)-  Given  that  Tsym  =  4  x  Tc  was  used  in  this  illustration,  the  spread  signal  PSD  oc¬ 
cupies  four  times  more  bandwidth  than  the  unspread  signal  when  measured  between  nulls. 
This  bandwidth  expansion  is  equivalent  to  what  is  commonly  called  the  DSSS  processing 
gain  which  is  given  by  Gp  =  Rc/Rsym  for  BPSK  data  modulation.  It  is  also  important 
to  note  that  the  peak  value  of  Pdsss  one-fourth  (l/Gp)  the  value  of  Pbpsk  (—6  dB  on  a 
decibel  scale).  Thus,  the  DSSS  signal  power  is  effectively  spread  across  a  wider  bandwidth 
than  the  original  unspread  signal. 

The  DSSS  receiver  essentially  despreads  and  estimates  communication  symbols  by 
repeating  transmitter  functions  in  reverse  order.  This  is  functionally  illustrated  in  the  block 
diagram  shown  in  Fig.  2.9.  The  received  DSSS  signal  enters  the  system  and  is  first  despread 
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Figure  2.8:  PSD  comparison  of  unspread  data  modulated  carrier  (dashed 
line)  and  final  spread  DSSS  waveform  (solid  line).  The  amount  of  spectral 
spread  is  dictated  by  processing  gain  GPp  =  Tgym/Tc  [25]. 


by  mixing  it  with  an  estimate  of  the  transmitted  spreading  waveform  c{t).  To  be  effective 
when  there  is  geographic  separation  between  the  transmitter  and  receiver,  the  receiver  must 
have  some  a  priori  knowledge  of  the  transmitted  spreading  code  and  thus  the  code  cannot 
be  purely  random.  The  despread  signal  is  bandpass  filtered  and  the  communication  symbols 
estimated  using  techniques  that  are  identical  to  systems  employing  no  spread  spectrum 
techniques. 


Figure  2.9:  Illustration  of  typical  DSSS  receiver  architecture  [25]. 
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During  DSSS  receiver  processing,  the  noise  and  other  interfering  signal  components 
can  be  thought  of  as  being  spread  versus  despread  given  they  do  not  contain  the  original 
spreading  modulation.  Thus,  the  composite  despread  waveform  consists  of  two  terms: 
1)  one  due  to  the  desired  transmitted  signal  which  returns  to  its  original  spectral  shape  and 
bandwidth,  and  2)  undesired  noise  and  interfering  signal  contributions  that  are  spectrally 
spread  with  reduced  peak  power  levels.  After  post-despread  filtering,  a  majority  of  the 
desired  signal  power  remains  (70%  to  90%  depending  of  the  filter  bandwidth)  while  only  a 
fraction  of  the  unwanted  noise  and  interfering  power  remain.  The  portion  of  desired  signal 
power  that  remains  and  the  percentage  of  undesired  power  that  is  rejected  is  a  function  of 
the  DSSS  processing  gain  Gp.  In  general,  the  SNR  improvement  due  to  the  despreading 
and  filtering  operation  is  proportional  to  Gp  [25] . 

For  a  DSSS  system  to  optimally  estimate  communication  symbols,  the  DSSS  receiver 
must  have  some  knowledge  of  the  received  signal  characteristics  across  the  signaling  do¬ 
mains:  time,  frequency,  space,  polarization,  and  code.  Depending  on  the  system,  some 
of  these  parameters  are  known  a  priori  while  others  must  be  estimated.  In  addition,  the 
estimated  parameters  may  change  over  time  and  must  be  constantly  tracked  and  updated. 
For  example,  a  communication  system  often  operates  at  a  pre-determined  center  frequency 
fc-  However,  frequency  variation  in  the  received  signal,  due  to  Doppler  frequency  shift 
fd  resulting  from  relative  transmitter-receiver  motion  and  imperfect  local  oscillator  behav¬ 
ior  in  both  the  transmitter  and  receiver,  dictates  that  the  receiver  employ  frequency  and/or 
phase  tracking  which  is  typically  accomplished  with  a  Phase-Lock  Loop  (PLL).  Next  to 
PLL  tracking  stability,  the  next  most  important  tracking  requirement  is  perhaps  the  abil¬ 
ity  to  reliably  generate  a  local  estimate  of  the  received  spreading  code.  The  next  section 
describes  one  common  tracking  method  used  for  DSSS  code  tracking. 

2.1.3  Delay -Lock  Loop  (DLL)  Code  Tracking.  Spreading  code  tracking  is  per¬ 
haps  the  most  important  aspect  in  a  DSSS  system.  While  the  transmitted  code  and  code 
parameters  are  generally  known  a  priori  by  the  receiver,  the  relative  time  offset  or  delay  of 
the  code  (sometimes  called  it  phase)  and  chip-to-chip  interval  variation  must  be  estimated 
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and  tracked  -  this  is  referred  to  as  code  tracking.  In  a  eommunieation  system,  a  laek  of  ef- 
feetive  eode  traeking  generally  results  in  ineffieient  despreading  of  the  reeeived  signal.  This 
results  in  a  lower  SNR  at  the  demodulator  input  and  poorer  bit  error  performanee  [15, 16]. 
In  applieations  other  than  eommunieations  that  also  employ  some  form  of  eode  traeking, 
e.g.,  spread  speetrum  based  navigation  or  radar  systems,  ineffeetive  eode  traeking  may 
result  in  poor  delay  estimates  whieh  ultimately  lead  to  poorer  position  estimates  [21]. 

To  traek  the  reeeived  PN  eoded  waveform,  the  DSSS  reeeiver  must  first  seareh  and 
find  the  relative  eode  position  using  a  proeess  ealled  aequisition.  This  proeess  provides 
an  initial  eourse  estimate  of  the  eode’s  temporal  position  and  is  used  as  a  starting  point 
for  more  aeeurate  traeking  by  the  eode  traeking  loop.  Onee  the  traeking  loop  reeeives  an 
initial  temporal  estimate  of  eode  position,  it  undergoes  a  proeess  ealled  pull-in,  where  it 
fine-tunes  the  estimated  eode  position.  Onee  it  has  adjusted  the  loop  enters  a  state  known 
a  traeking,  where  it  is  providing  a  suffieiently  aeeurate  eode  estimate,  in  terms  of  temporal 
position,  to  effeetively  despread  the  reeeived  DSSS  waveform.  The  required  eode  traeking 
aeeuraey  varies  with  system  requirements. 

Code  traeking  is  generally  performed  using  PLL  teehniques.  The  most  eommonly 
used  arehiteeture  employs  two  parallel  branehes,  ineluding  one  representing  an  early  (ad- 
vaneed)  version  of  the  eurrent  eode  estimate  and  the  other  representing  a  late  (delayed) 
version  of  the  eurrent  eode  estimate.  This  arehiteeture  is  known  as  a  Delay-Look  Loop 
(DLL)  [32]  and  is  the  foous  of  this  researoh.  DLL  implementations  ean  be  eategorized 
as  being  either  eoherent  or  non-eoherent.  A  ooherent  DLL  uses  knowledge  of  the  reeeived 
signal  oarrier  frequenoy  and  phase  to  perform  traeking.  Therefore,  the  signal  may  be  down- 
eonverted  with  the  eode  traeking  ooourring  at  baseband.  However,  in  many  situations  it  ean 
be  diffioult  to  estimate  and  traek  the  oarrier  phase  without  first  traeking  the  eode  for  ef- 
feotive  despreading.  For  this  reason,  the  analysis  in  this  dissertation  oonoentrated  on  the 
non-eoherent  DLL,  whieh  typieally  operates  at  an  intermediated  frequenoy  (IF)  and  does 
not  require  or  assume  knowledge  of  the  earrier  phase. 
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Figure  2.10:  Typical  elements  of  non-coherent  Delay -Lock  Loop  (DLL)  [25]. 

The  typical  elements  for  implementing  a  non-coherent  DLL  are  shown  in  Fig.  2.10. 
The  two  distinct  parallel  branches  following  the  RF  BPF  are  referred  to  as  the  early  and 
late  gate  branches.  These  two  branches  are  functionally  equivalent  beginning  with  the  left¬ 
most  despreading  mixer  and  ending  with  the  right-most  LPF.  What  distinguishes  the  two 
branches  is  that  the  estimated  code  from  the  code  generator  is  either  delayed  (arrives  late) 
or  advanced  (arrives  early)  by  some  fixed  amount  before  being  applied  to  the  despreading 
mixer.  One  intuitive  way  to  understand  DLL  functionality  is  to  observe  how  the  time  and 
frequency  domain  responses  of  a  given  input  signal  are  modified  as  the  signal  progresses 
through  the  DLL.  Given  the  early  and  late  gate  branches  are  functionally  equivalent  except 
for  a  code  timing  offset,  the  process  is  illustrated  using  Cg/z  (t)  to  represent  either  the  early 
or  late  gate  code. 

Received  signal  Srit)  =  dit)cit)cos  [27r(/c  -f  fd)t],  where  fc  is  the  transmitted  car¬ 
rier  frequency  and  fd  is  Doppler  frequency  shift,  is  first  passed  through  the  DLL  RF  BPF. 
Given  an  RF  BPF  bandwidth  of  Wrf  =  ‘^RcUp  and  an  arbitrary  portion  of  Srit)  that  spans 
two  chip  intervals,  the  DLL  RF  BPF  produces  the  filtered  time  and  frequency  domain  re¬ 
sponses,  Slit)  and  Siif),  shown  in  Fig.  2.11  and  Fig.  2.12,  respectively.  Selection  of  the 
RF  BPF  bandwidth  is  usually  related  to  the  main  spectral  response  of  the  received  signal. 
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This  allows  the  filter  to  pass  most  of  the  desired  signal  energy  while  blocking  undesired 
signal  and  noise  responses  falling  outside  the  band  of  interest.  Notice  in  Fig.  2.1 1  that  less 
than  ideal  filtering  removes  the  ideal  rectangular  shape  across  the  chip  boundary  and  that 
only  the  main  spectral  response  remains  in  Fig.  2.12. 

The  remaining  signal  is  then  multiplied  by  an  offset  version  of  the  baseband  PN  code. 
The  time  domain  result  S2{t)  is  shown  in  Fig.  2.13  and  the  frequency  domain  result  S2{f) 
is  in  Fig.  2.14.  The  important  observation  at  this  point  occurs  in  the  frequency  domain, 
where  one  can  see  that  spectral  lines  are  created  at  harmonics  of  the  code  repetition  rate, 
Rc.  The  magnitude  of  the  spectral  line  a  is  used  for  code  tracking.  If  the  code  is  perfectly 
tracked,  this  spectral  line  will  have  the  same  magnitude  in  both  the  early  and  late  branches 
of  the  DLL.  If  it  is  not  perfectly  tracked,  then  the  magnitude  of  this  spectral  line  in  the  two 
branches  will  differ. 

To  recover  the  desired  spectral  line  at  fc,  the  signal  is  passed  through  the  IF  BPF 
The  bandwidth  of  this  filter  should  be  as  narrow  as  possible  to  remove  unwanted  noise  and 
harmonic  effects,  while  being  wide  enough  to  ensure  sufficient  desired  signal  energy  passes 
given  that  received  signal  frequency  (fc  +  fd)  is  not  precisely  known.  The  resultant  time 
and  frequency  domain  responses  of  the  DSSS  signal  at  the  DLL  IF  BPF  output,  Ss(t)  and 
Ss(f)  are  shown  in  Fig.  2.15  and  Fig.  2.16,  respectively. 

The  magnitude  (envelope)  of  the  DLL  IF  BPF  output  effectively  represents  a  measure 
of  correlation  between  the  received  PN  coded  waveform  and  the  DLL  estimate  of  the  same. 
The  envelope  is  extracted  using  a  basic  energy  detection  process  comprised  of  a  squaring 
operation  followed  by  lowpass  filtering  (LPF).  The  result  of  applying  this  operation  to 
the  signal  in  Fig.  2.15  yields  the  time  and  frequency  domain  responses  shown  in  Fig.  2.17 
and  Fig.  2.18,  respectively. 

The  filtered  difference  between  the  early  and  late  gate  energy  detector  outputs,  or 
discriminator  output,  provides  a  measure  of  relative  code  offset  between  the  received  and 
internally  generated  codes.  The  loop  filter  design  is  highly  application  specific  and  aims 
to  maximize  overall  DLL  stability.  The  discriminator  output  controls  the  numerically  con- 
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Figure  2.11:  Time  domain  response  of  DSSS  signal  at  the  DLL  RF  BPF 
output. 


fc  —  Rc  fc  fc  +  Rc 

Frequeney 


Figure  2.12:  Frequeney  domain  response  of  DSSS  signal  at  the  DLL  RF 
BPF  output. 
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Figure  2.13:  Unfiltered  time  domain  response  of  despread  DSSS  signal 
after  multiplieation  by  c{t)  under  perfeet  eode  traeking  eonditions. 


/c  ~  Rc  fc  fc  +  Rc 

Frequeney 


Figure  2.14:  Unfiltered  frequeney  domain  response  of  despread  DSSS 
signal  after  multiplieation  by  c{t)  under  perfeet  eode  traeking  eonditions. 
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Figure  2.15:  Filtered  time  domain  response  of  despread  DSSS  signal  at 
DLL  IF  BPF  output  under  perfeet  eode  traeking  eonditions. 


Frequeney 


Figure  2.16:  Filtered  frequeney  domain  response  of  despread  DSSS  sig¬ 
nal  at  DLL  IF  BPF  output  under  perfeet  eode  traeking  eonditions. 
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trolled  oscillator  (NCO)  which  drives  the  PN  code  generator.  If  the  current  state  of  the 
internal  code  generator  perfectly  matches  the  received  code,  the  discriminator  output  is 
zero  and  no  change  occurs  in  the  NCO.  Whenever  the  internal  code  generator  state  does 
not  perfectly  match  the  received  code,  a  non-zero  signed  discriminator  output  is  produced 
and  the  DLL  code  generator  advances  or  delays  its  output  accordingly  (changes  the  code 
phase).  For  example,  if  the  early  gate  correlates  more  closely  with  the  received  signal  the 
discriminator  output  becomes  positive  and  the  NCO  frequency  increases.  Likewise,  if  the 
late  gate  correlates  more  closely  with  the  received  signal  the  discriminator  output  becomes 
negative  and  the  NCO  frequency  decreases.  Thus,  the  DLL  is  constantly  adapting  to  signal 
and  channel  conditions.  This  adaptation  is  desirable  when  compensating  for  actual  varia¬ 
tions  in  Sr{t)  such  as  changes  in  carrier  frequency  fc  and  Doppler  frequency  fd-  However, 
DLL  tracking  variation  due  to  noise  and/or  interfering  signals  is  generally  undesirable.  The 
DLL  output  discriminator  response  for  changing  signal  conditions  is  generally  described  as 
an  S-curve.  An  ideal  S-curve  response  is  shown  in  Fig.  2.19  for  the  case  where  no  channel 
noise  or  interfering  signals  are  present. 
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Time 


Figure  2.17:  Time  domain  response  of  DSSS  signal  at  the  energy  detee- 
tor  output  showing  that  double  frequeney  and  higher-order  harmonies  have 
been  suppressed. 


Frequeney 


Figure  2.18:  Frequeney  domain  response  of  DSSS  signal  at  the  energy 
deteetor  output  showing  that  double  frequeney  and  higher-order  harmonies 
have  been  suppressed. 
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Discriminator  Output 


Figure  2.19:  Ideal  S-eurve  response:  DLL  diseriminator  output  versus 
eode  phase  offset  between  the  reeeived  and  DLL  estimated  eodes  [25]. 
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2.2  Optimization  Techniques 

Optimization  techniques  attempt  to  maximize  or  minimize  the  response  of  a  math¬ 
ematical  model  that  represents  a  given  physical  system.  The  mathematical  function  to  be 
maximized  or  minimized  is  called  the  objective  function,  which  represents  the  system  re¬ 
sponse  to  a  particular  set  of  decision  variables.  The  relationship  between  the  objective 
function  and  the  decision  variables  is  described  by  model  parameters  and  constraints.  In 
general,  there  may  be  many  combinations  of  decision  variables  that  satisfy  the  constraints, 
the  collection  of  which  is  known  as  the  feasibility  region.  Therefore,  an  optimized  solution 
for  the  system  is  the  “best  choice”  of  decision  variables  that  fall  within  the  feasibility  region 
-  the  particular  definition  of  “best  choice”  is  problem- specific  [18]. 

2.2.1  Genetic  Algorithm  (GA).  There  are  many  optimization  techniques  which 
may  be  applied  to  a  particular  problem,  with  each  having  its  own  strengths  and  weaknesses. 
One  large  class  of  techniques  is  known  as  heuristic  methods  which  generally  attempt  to  find 
a  “good”  solution  without  necessarily  striving  to  guarantee  optimality  [18].  While  heuristic 
techniques  tend  to  be  very  problem- specific,  there  are  a  handful  of  metaheuristic  methods 
which  may  be  applied  to  a  broader  range  of  problems.  The  GA  is  one  popular  class  of 
metaheuristics  which  simulate  the  biological  evolution  process  by  describing  the  decision 
parameters  using  a  binary  string  called  a  gene.  The  GA  process  essentially  mates  genes 
from  an  available  population  and  retains  the  strongest  offspring  for  subsequent  mating  in 
the  new  population.  Some  parents  are  retained  in  the  new  population  and  the  possibility  of 
mutation  within  offspring  permitted  [18].  The  GA  process  can  be  summarized  as  follows: 

1 .  The  initial  GA  population  is  randomly  generated  and  genes  formed  using  the  process 
illustrated  in  Fig.  2.20.  The  fitness  of  each  generated  gene  is  then  calculated  to  ensure 
it  is  in  the  feasibility  region.  If  in  the  feasibility  region  it  is  retained,  else,  it  is  rejected 
and  another  gene  replaces  it. 

2.  Some  number  of  most  fit  (m)  and  least  fit  (/)  genes,  for  (I  +  m)  a  multiple  of  two, 
are  selected  from  the  population  to  serve  as  parents.  The  selected  parents  are  then 
randomly  paired  for  mating  to  create  {m  +  1)/ 2  parent  pairs. 
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Figure  2.20:  GA  genes  assembled  for  the  initial  population  using  two 
parameters  with  randomly  generated  binary  digit  values. 


3.  The  binary  digits  of  the  parent  genes  are  eompared  to  generate  offspring  as  shown 
in  Fig.  2.21.  When  binary  digits  in  the  parents  mateh,  their  eorresponding  values  are 
passed  direetly  to  the  offspring.  When  the  binary  digits  differ  (indieated  by  the  x 
elements  in  the  figure),  the  offspring  digit  values  are  randomly  assigned  with  equal 
probability. 

4.  In  addition,  eaeh  digit  in  the  seleeted  offspring  is  subjeeted  to  some  likelihood  of 
mutation,  i.e.,  there  is  some  ehanee  that  its  final  value  is  eomplemented  as  illustrated 
in  Fig.  2.22. 

5.  A  new  gene  population  is  formed  by  retaining  some  of  the  previously  mated  parents 
and  their  offspring. 

6.  The  seleetion-mating-seleetion  proeess  in  Step  2  through  Step  5  is  repeated  for  sev¬ 
eral  iterations  (generations)  until  an  exit  eriterion  is  met. 
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Figure  2.21:  GA  offspring  generation  process  using  a  chosen  parent  gene 
pair  from  the  population.  The  boxed  columns  denote  dissimilar  parent  val¬ 
ues  where  offspring  values  are  randomly  assigned.  Matching  parent  values 
in  unboxed  columns  are  assigned  directly  to  the  offspring. 
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Figure  2.22:  GA  offspring  mutation  process.  Each  offspring  digit  is  sub¬ 
jected  to  some  likelihood  of  mutation. 
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2.2.2  Response  Surface  Methodology  (RSM).  RSM  is  a  statistically -based  al¬ 
ternative  to  GA  optimization  and  proven  itself  in  industrial  experimentation  and  typieally 
involves  quality  eontrol.  However,  it  is  believed  that  this  same  approaeh  may  be  useful  in 
other  types  of  engineering  applieations,  ineluding  waveform  design.  Given  it  has  a  solid 
foundation  in  linear  systems  theory,  it  is  readily  aeeessible  for  oommunieations  design  and 
eonsists  of  the  following  eomponents  [22] : 

1 .  Design  of  Experiments  (DOE),  Seetion  2.2.2. 1 :  A  proeess  for  setting  up  experiments 
to  ensure  proper  eolleetion  of  data.  If  designed  eorreetly,  the  experiment  provides 
data  that  may  be  used  for  linear  regression  modeling  and/or  Analysis  of  Varianee 
(ANOVA)  proeessing  [7,22]. 

2.  Einear  Regression  Modeling,  Seetion  2. 2. 2.2:  Used  in  RSM  to  estimate  the  shape  of 
a  response  surfaee.  If  the  response  surfaee  eontains  eurvature  over  the  tested  region, 
the  surfaee  model  may  be  used  to  solve  for  an  optimal  point.  If  the  surfaee  does  not 
eontain  eurvature,  the  regression  model  may  still  be  used  to  seareh  for  a  region  that 
is  quadratie  [22] . 

3.  Analysis  of  Varianee  (ANOVA),  Seetion  2. 2. 2. 3:  A  statistieal  analysis  tool  used  to 
statistieally  eompare  the  mean  values  of  eolleeted  data.  The  ANOVA  proeess  is 
elosely  related  to  linear  regression  modeling  and  may  be  used  to  validate  a  given 
regression  model,  or  to  determine  the  signifieanee  of  various  model  elements,  in¬ 
eluding  the  presenee  of  eurvature  in  the  response  surfaee.  Conelusions  drawn  from 
the  ANOVA  proeess  are  based  on  the  assumption  that  the  residual  data  is  normally 
distributed  and  that  the  population  varianees  are  equal  [7,22]. 

4.  Residual  Testing,  Seetion  2. 2. 2.4  through  Seetion  2. 2. 2.7:  Provides  methods  for  test¬ 
ing  ANOVA  residual  data  distributions  to  verify  whether  or  not  the  normality  as¬ 
sumptions  are  satisfied.  If  the  assumptions  are  not  satisfied,  various  transformation 
teehniques  ean  be  applied  to  remedy  violations  [14, 22]. 

5.  Comparative  Testing,  Seetion  2. 2. 2. 8:  The  ANOVA  null  hypothesis  is  that  the  mean 
values  of  the  tested  populations  are  equal.  If  the  null  hypothesis  is  rejeeted  (means 
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are  not  equal),  at  least  one  of  the  population  means  is  statistieally  different.  In  and 
of  itself,  the  ANOVA  proeess  does  not  provide  insight  as  to  which  population(s)  are 
different.  Thus,  additional  comparative  testing  is  required  to  make  this  determina¬ 
tion  [14,22]. 

2.2.2. 1  Design  of  Experiments  (DOE).  DOE  is  a  process  that  allows  exper¬ 
imenters  to  craft  trials  that  can  statistically  characterize  or  optimize  a  process  or  a  system. 
The  system  under  test  may  often  be  modeled  as  a  linear  process  given  by  [22], 

Y  =  X/3  +  e  ,  (2.8) 

where  Y  is  the  system  response,  X  is  a  matrix  representing  combinations  of  input  parame¬ 
ters,  (3  is  the  system  model,  and  e  is  random  error.  Several  goals  may  be  accomplished  with 
a  design  such  as  this.  First,  one  may  determine  if  the  system  model  (3  adequately  describes 
true  system  behavior.  Second,  one  may  be  able  to  establish  a  model  estimate  (3.  Third,  if 
P  can  be  determined,  one  can  gain  knowledge  to  optimize  response  Y  given  parametric 
variation  in  X  [22] . 

The  input  parameters  in  matrix  X  are  often  expressed  in  terms  of  coded  variables  [22] . 
This  is  done  by  mapping  between  natural  variables  and  coded  parameters  in  X.  The  input 
variables  in  this  dissertation  will  generally  be  the  SMSE  waveform  design  variables  Nf 
and  A/.  Given  these  variables,  a  representative  mapping  from  SMSE  parameters  (natural 
units)  to  DOE  coded  units  is  shown  in  Table  2.1. 

Table  2.1:  Representative  mapping  from  SMSE  variables 
(Nf,  A/)  to  DOE  coded  variables  (xi,  X2). 


Factor  Eevel 

Nf 

Xi 

A/ 

X2 

Eow 

8 

-1 

17 

-1 

High 

32 

1 

15 

1 

32 


By  way  of  illustrating  the  use  of  (2.8)  with  the  mapping  in  Table  2.1,  eonsider  a  first- 
order  model  with  interaction.  In  this  case,  each  element  in  the  vector  Y  can  be  expressed 
as 


1/21  =  /3o  +  /5lXi  +  (32X2  +  (3i2XiX2  +  6  .  (2.9) 

Using  this  model  to  test  the  condition  Nf  =  32  and  A/  =  17,  the  corresponding  mapped 
values  of  xi  =  1  and  0:2  =  —1  from  Table  2.1  are  used  such  that  (2.9)  is  expressed  as 

//  =  /?o  +  A(l)  +  /^2(— 1)  +  /3i2(1)(— 1)  +  e. 

A  design  that  would  test  all  possible  combinations  of  (iV/,  A/)  as  given  in  Table  2.1  is 
called  a.  full  factorial  design  [22]  and  would  include  all  elements  Y  given  by 


Uii  =  /?o  + /3i(~l)  + /32(~1)  + /3i2(~1)(~1)  +  e 
yi2  =  (3o  +  (3i{—l)  +  (32{f)  +  (3i2{—f){l)  +  e 
1/21  =  (3o  +  (diil)  +  (32^—f)  +  (3i2{X){—f)  +  e 
1/22  =  /^o  + /3i(l)  + /52(1)  + /3i2(1)(1)  +  e  • 


This  set  of  equations  may  be  expressed  in  the  matrix  form  given  by  (2.8)  using 


1 


and 


X 


FF  — 


1 

1 

1 

1 


-1 

-1 

1 

1 


-1 

1 

-1 

1 


-1 

-1 

1 


(2.11) 


/^o 

A 

A 

Pl2 


(2.12) 


where  subseript  values  in  Y  and  jS  eorrespond  to  ordered  eoded  variables  (xi,  X2)  and  Xff 
denotes  a  full  faetorial  design. 

If  experimental  replieation  is  desired  for  statistieal  analysis,  which  is  often  the  case 
for  optimization,  the  input  parameters  in  XpF  may  be  repeated.  This  is  accounted  for  using 


X  = 


Xff 

Xff 


X 


FF 


(2.13) 


where  X  now  represents  the  total  experiment  and  Xff  represents  one  repetition  of  the  full 
factorial  design  [22].  As  provided,  the  expressions  in  (2.8),  (2.10),  (2.11),  and  (2.12)  pro¬ 
vide  a  full  experimental  description  that  may  be  used  directly  to  analyze  system  response 
Y  [22]. 


2.2.22  Linear  Regression  Modeling.  Once  an  experiment  has  been  com¬ 
pleted  according  to  (2.8),  the  system  response  Y  (collected  data)  may  be  used  to  generate  a 
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model  estimate  (3.  This  may  be  done  using  a  linear  regression  proeess  by  solving  the  least 
squares  normal  equation  given  by  [22] 

=  (X^X)  X'^Y  ,  (2.14) 

with  the  predieted  values  of  Y  then  given  by 

Y  =  X/3  .  (2.15) 

The  resultant  residuals  e  (error)  are  then  ealeulated  as  the  differenee  between  the  eollected 
and  predieted  values  of  Y  aecording  to  [22] 

e  =  Y-Y.  (2.16) 

If  the  model  estimate  P  provides  a  good  predietion  of  response  Y  given  input  X,  then 
the  residual  elements  in  e  will  be  normally  distributed  random  variables.  Furthermore,  the 
varianee  of  e  under  eaeh  eondition  in  X  should  be  approximately  equal  [22].  Tests  for 
quantifying  normality  and  varianee  equality  among  the  residual  elements  are  deseribed  in 
Section  2.2.2. 5  and  Section  2. 2. 2. 6,  respectively.  If  the  model  estimate  /3  indicates  there 
is  curvature  in  the  response  surface  Y,  it  may  be  used  to  optimize  Y,  i.e.,  $  may  be  used 
to  specify  optimal  input  parameters  in  X  to  either  maximize  or  minimize  response  surface 
Y  [22]. 

Two  /3  models  are  particularly  useful  for  RSM:  the  second-order  model  and  the  first- 
order  model  [22] .  The  second-order  model  is  given  by 


y  —  Po  3-  PiXi  +  P2X2  +  (3i2XiX2  +  (3iix\  -\-  P22X2  ■  (2.17) 

If  the  estimated  model  parameters  from  (2.17)  adequately  describes  the  data  (as  determined 
by  the  ANOVA),  then  the  surface  is  assumed  to  contain  curvature  and  optimal  input  param- 
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eters  xi  and  X2  may  be  calculated.  Using  elements  expressed  by  (2.17),  the  second-order 
model  may  be  expressed  in  matrix  form  as 

Y  = /3o  +  X'^b  +  X'^BX  ,  (2.18) 


where 


X2 


J2  _ 

and 

A  stationary  point  for  the  second-order  model  is  given  by  [22] 

Xs  =  -^B-ib.  (2.19) 

The  derivation  of  stationary  point  Xg  is  accomplished  by  setting  the  derivative  of  the  es¬ 
timated  model  P  to  zero  and  solving  for  X  [22].  Therefore,  the  stationary  point  must  be 
tested  to  determine  if  it  creates  a  maximum,  minimum,  or  saddle  point  in  response  surface 
Y.  This  is  done  using  Eigenvalue  analysis.  If  all  Eigenvalues  of  B  are  positive,  Xg  mini- 
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mizes  Y.  If  all  Eigenvalues  are  negative,  Xg  maximizes  Y.  If  the  Eigenvalues  signs  for  B 
differ,  the  point  Xg  eorresponds  to  a  saddle  point  on  Y  [22] . 

If,  however,  when  using  the  seeond-order  model  of  (2.17)  the  resultant  ANOVA  does 
not  determine  the  response  surfaee  Y  eontains  eurvature,  a  first-order  model  of  the  follow¬ 
ing  form  may  be  used  [22] 


2/  —  /?o  +  PiXi  +  (^2X2  ■ 

In  this  ease,  the  path  of  steepest  aseent/deseent,  given  by 


(2.20) 


Axi  Ax2 
Pi  P2 


(2.21) 


may  be  followed  in  subsequent  experiments  until  the  seeond-order  model  beeomes  appro¬ 
priate.  Then,  a  stationary  point  may  be  found  by  (2.19)  [22]. 


2.2.23  Analysis  of  Variance  (ANOVA).  The  ANOVA  proeess  is  the  heart 
of  DOE  optimization.  ANOVA  is  a  statistieal  teehnique  whieh  uses  the  sample  varianees 
of  a  data  set  to  test  the  impaet  of  input  parameters  on  an  output  response  variable.  Speeif- 
ieally,  given  an  experimental  design  X  and  model  P,  the  ANOVA  proeess  uses  the  sample 
varianees  to  test  the  null  hypothesis,  i.e.,  “Are  the  means  of  different  eonditions  deseribed 
by  the  model  equal?”  Thus,  rejeetion  of  the  null  hypothesis  implies  that  at  least  one  of  the 
means  are  different  [22]. 

There  are  many  ways  to  interpret  a  null  hypothesis  rejeetion.  In  this  analysis,  the 
ANOVA  is  primarily  used  in  two  ways.  Eirst,  rejeeting  the  null  hypothesis  with  respeet  to 
the  entire  regression  model  signifies  that  P  adequately  deseribes  the  proeess  [22].  Seeond, 
the  ANOVA  is  used  to  evaluate  the  individual  input  parameters  of  X.  In  this  ease,  rejeeting 
the  null  hypothesis  for  an  given  input  parameter  indieates  that  parameter  has  a  signifieant 
effeet  on  response  surfaee  Y  and  should  be  ineluded  in  the  model  for  optimization  [22] .  As 
a  result,  ANOVA  provides  insight  into  the  response  surfaee  shape  by  ineluding  appropriate 
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model  terms  while  eliminating  those  that  are  inappropriate.  The  ANOVA  proeess  begins  by 
assuming  a  model  for  the  test  data  [7].  For  this  illustration,  the  assumed  model  is  second- 
order  and  expressed  as  [22] 


Uijk  —  /^O  +  +  (^2X2^  +  (3l2XlfX2j 

+  f^nxl.  +  (d22x\.  +  Cijk  1  (2.22) 

where  y  is  the  response  variable,  (xi,  X2)  are  coded  input  variables,  the  (3  terms  represent 
regression  coefficients  and  e  represents  the  error.  The  i  and  j  subscripts  correspond  to  index 
values  of  xi  and  X2,  respectively.  The  k  subscript  represents  the  experimental  run  number. 
The  expression  in  (2.22)  may  be  expressed  more  compactly  in  matrix  form  as 

Y  =  X/3  +  e  .  (2.23) 

Essentially,  the  ANOVA  process  performs  a  least-squares  fit  of  the  data  to  the  model  by 
applying  [7,22] 


(5  =  (X^X)  ^  X^Y  .  (2.24) 

Substituting  this  estimate  for  (3  into  (2.23)  yields  a  regression  model  of  the  form 

Y  =  X/3  =  X  (X^X)  X^Y  ,  (2.25) 

with  the  resultant  error  given  by 


e  =  Y  -  Y  .  (2.26) 

After  the  least-squares  fit  to  the  data,  the  sample  variance  is  partitioned  into  sub¬ 
spaces  corresponding  to  the  main  effects,  interactions  and  error.  The  main  effects  variance 
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and  interaction  variances  are  compared  to  the  error  variance  as  part  of  a  significance  test 
to  determine  their  relative  importance  to  the  model.  More  specifically,  assuming  that  the 
underlying  model  error  is  Gaussian,  the  ratio  of  a  given  factor’s  sample  variance  to  the 
error’s  sample  variance  will  be  distributed  according  to  the  F  distribution.  Tabulated  F 
distribution  values  may  then  be  used  to  provide  probabilities  that  the  given  factor  under 
test  is  significant  and  should  be  included  in  the  model.  The  larger  the  value  of  test  statis¬ 
tic  Fq,  the  more  reasonable  it  becomes  to  reject  the  null  hypothesis  of  equal  means.  An 
alternative  approach  uses  a  p- value  which  represents  the  probability  that  the  statistic  Fq  is 
the  result  of  a  population  that  corresponds  to  the  null  hypothesis.  If  the  p- value  is  small,  it 
casts  doubt  on  the  null  hypothesis.  Therefore,  ANOVA  provides  a  tool  for  deciding  which 
factors  and  interactions  are  most  important  in  a  given  system  model,  a  task  which  is  often 
left  to  intuition  [7,22].  More  detailed  information  regarding  ANOVA  processing  can  be 
found  in  [22]. 

2. 2. 2.4  Normality  and  Variance.  The  normality  and  equal  variance  as¬ 
sumptions  that  were  noted  in  Section  2. 2. 2. 2  and  Section  2. 2. 2. 3  are  important  to  both  lin¬ 
ear  regression  modeling  and  ANOVA  processing  [22].  The  following  subsections  address 
each  of  these  conditions. 

2.2.2. 5  Normality  Testing.  The  residuals  e  are  assumed  to  be  normally 
distributed  [22]  in  linear  regressing  modeling  and  ANOVA  processing.  A  visual  test  such 
as  a  probability  plot,  or  a  numerical  method  such  as  the  Shapiro- Wilks  test,  can  be  used  for 
testing  normality  [14]. 

A  representative  normal  probability  plot  of  residuals  is  shown  in  Fig.  2.23.  This  is 
essentially  a  plot  of  the  residuals  (x-axis)  versus  the  cumulative  probability  of  the  normal 
distribution  (p-axis).  If  the  residuals  appear  to  lie  on  a  straight  line  they  are  distributed 
approximately  normally  [14,22]. 

The  Shapiro-Wilks  test  provides  a  significance  test  for  normality  [14].  In  this  case, 
the  null  hypothesis  is  that  the  data  is  normally  distributed.  To  execute  the  test,  the  residuals 
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Residuals 

Figure  2.23:  Normality  test  using  residuals  versus  eumulative  probability 
of  normal  distribution.  If  residuals  lie  along  a  straight  line  the  data  may  be 
considered  normally  distributed  [22] . 

e  are  ordered  such  that  ei  <  62  <  ■  ■  •  <  where  Ne  represents  the  total  number  of 
residuals.  Test  statistic  W  is  then  calculated  using 


i=l 

Ne  even 
Ne  odd  , 

k 

b  =  ^  CLN^-i+l  (eAT^-i+l  -  Gi)  , 
i=l 


(2.27) 


(2.28) 


(2.29) 


,  (2.30) 

where  e  in  (2.27)  is  the  residual  mean  and  tabulated  values  of  in  (2.29)  are  provided 

in  [14].  The  resultant  value  of  W  from  (2.30)  is  then  mapped  to  a  tabulated  p-value  which 
can  be  found  in  [14].  Given  that  the  null  hypothesis  is  normally  distributed  data,  a  small  p- 
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value  from  the  Shapiro-Wilks  easts  doubt  on  this  assumption  and  causes  the  null  hypothesis 
to  be  rejected  [14]. 

2. 2. 2. 6  Variance  Testing.  Variance  equality  testing  can  be  done  using  a 
Bartlett’s  significant  test.  In  this  case,  the  null  hypothesis  is  that  the  residual  variances 
under  all  experimental  conditions  are  equal. 

Assuming  there  are  a  variances  to  compare  with  each  having  n*  samples,  the  total 
number  of  samples  N^,  is  given  by 


Ne  =  '^ni.  (2.31) 

i=l 

To  test  the  null  hypothesis  that  af  =  =  ■■■  =  al,  the  test  statistic  Xo  is  calculated  by 


S'!  = 


Nf,  —  a 


(2.32) 


i=l 


q=  {Ne-  a)  logio  “  1)  logic  i^i)  ,  (2.33) 

i=l 


C  =  1  + 


1 

3(a-  1) 


(Ne  -  !)■' 


(2.34) 


xl  =  2.3026  (2.35) 

where  Sf  represents  the  sample  variance  of  the  condition  in  the  experiment  [22].  The 
test  statistic  Xo  is  used  to  find  a  tabulated  p-value  from  a  xi-i  distribution.  If  the  p-value 
obtained  from  a  xl-i  distribution  is  small,  this  suggests  that  the  null  hypothesis  should  be 
rejected  and  the  data  does  not  meet  the  variance  equality  condition  [22]. 

2. 2. 2. 7  Statistical  Transformation.  If  the  normality  conditions  and/or  vari¬ 
ance  equality  conditions  in  Section  2. 2. 2. 5  and  Section  2. 2. 2. 6  are  not  satisfied,  it  may  be 
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possible  to  transform  the  data  so  the  eonditions  are  met.  The  Box-Cox  transformation  is 
one  eommon  transformation  for  aeeomplishing  this.  The  Box-Cox  transformation  is  given 
by  [22] 


T 


A 

In  (Y) 


A  7^  0 
A  =  0  , 


(2.36) 


where  A  G  i?.  As  indieated,  data  Y  may  be  transformed  using  a  range  of  A  values.  If  a 
value  of  A  exists  sueh  that  the  residuals  of  the  transformed  data  T  meet  the  requirements 
for  normality  and  equality  of  varianee,  then  linear  regression  and/or  the  ANOVA  may  be 
performed  on  the  transformed  data,  T  [22] .  However,  if  the  residuals  from  T  still  does  not 
satisfy  normality  and  varianee  equality  eonditions,  the  ranks  of  the  data  ean  be  analyzed  as 
part  of  the  ANOVA  proeess.  When  ranks  are  used,  the  test  is  known  as  the  Kruskal-Wallis 
test  [22]. 


2. 2.2.8  Comparative  Testing.  The  ANOVA  proeess  is  a  hypothesis  test 
that  provides  information  about  the  statistieal  data  means  for  eonditions  under  test.  The 
null  hypothesis  is  that  all  means  are  equal  and  its  rejeetion  indieates  that  at  least  one  mean 
differs  from  the  others.  However,  in  and  of  itself  the  ANOVA  proeess  does  not  provide 
an  indieation  of  whieh  mean(s)  is  different.  Comparative  testing  ean  be  performed  to  de¬ 
termine  this.  In  addition  to  providing  information  about  the  ANOVA  results,  eomparative 
testing  is  also  important  to  quantify  results  of  the  optimization  proeess  [22] . 

One  visual  tool  for  eomparing  population  statisties  for  random  variables  is  the  box 
and  whisker  plot  as  illustrated  in  Fig.  2.24  [22].  This  plot  shows  statistieal  properties  of 
the  eorrelation  degradation  metrie  Coeg  under  various  SMSE  parameter  eombinations.  The 
box  and  whisker  plot  is  interpreted  as  follows  for  a  given  SMSE  parameter  eombination: 
2)  the  box  midline  represents  the  median  value,  2)  the  top  and  bottom  box  edges  represent 
quartiles  for  the  25*^  and  75*^  percentiles  of  the  populations,  and  3)  the  extreme  “whisker” 
ends  represent  the  minimum  and  maximum  population  values.  The  box  and  whisker  repre¬ 
sentation  reveals  general  trends  about  the  data  means  and  data  spread  [22]. 
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Figure  2.24:  Representative  box  and  whisker  plot  showing  population 
medians  (box  midline),  quartiles  for  25*^  and  75*^  percentiles  (box  top  and 
bottom),  and  population  extreme  values  (whisker  ends)  [22]. 

To  formally  compare  the  means  of  two  normally  distributed,  equal  variance  popula¬ 
tions,  a  f-test  may  be  used  according  to  the  following  [22].  Given  there  are  n*  samples  in 
the  population,  with  being  the  mean  value  for  the  the  population,  test  statistic  to 
for  the  two  mean  case  can  be  calculated  as  follows: 


123456789  10 

SMSE  Parameter  Combination 


fii  1 


(2.37) 


k=l 


2  (ni  -  1)  Si  +  (ria  -  1)  Si 


s;  = 


77,1+77,2-2 


(2.38) 


to  =  .  (2.39) 

Sp^/  —  +  — 

Py  ni  712 

Test  statistic  to  from  (2.39)  is  used  to  find  a  tabulated  p- value  from  a  fni-i-n2-2  distri¬ 
bution.  If  the  p- value  obtained  from  the  tn-^+n2-2  distribution  is  small,  this  suggests  that  the 
null  hypothesis  should  be  rejected  and  the  means  and  ^2  do  indeed  differ  [22]. 
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To  compare  means  for  more  than  two  populations,  the  Least  Signifieant  Differenee 
(LSD)  test  may  be  used.  Assuming  there  are  a  total  means  to  eompare  with  eaeh  based  on 
Ui  samples,  the  total  number  of  samples  Ng,  is  given  by 

a 

Ne  =  '^ni  .  (2.40) 

i=l 

In  this  ease,  the  LSD  test  for  a  full-faetorial  design  is  given  by 


LSD  =  t 


2  e'^e 

2Ae-a(/ 


(2.41) 


where  a  is  the  signifieanee  level  and  u  =  —  a  is  the  number  of  degrees  of  freedom  [22] . 

When  eomparing  any  two  populations  in  the  experiment,  the  null  hypothesis  of  equal  means 
(Vi  =  Vj)  would  be  rejeeted  if  the  means  differ  by  more  than  the  LSD  [22]. 
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III.  Methodology 

In  many  respects,  the  goal  of  this  dissertation  is  to  develop  and  demonstrate  a  set  of  tools 
that  waveform  design  engineers  can  use  to  design  efficient  communication  waveforms. 
One  of  the  basic  tools  adopted  for  this  research  is  the  Spectrally  Modulated,  Spectrally 
Encoded  (SMSE)  framework  as  described  in  Section  2.1.1.  The  SMSE  framework  has 
been  successfully  used  to  mathematically  describe  various  types  of  Orthogonal  Erequency 
Division  (OEDM)  waveforms.  The  quantitative  nature  of  the  SMSE  framework  makes 
it  well-suited  for  applying  design  techniques  from  disciplines  outside  of  communications 
and/or  signal  processing.  More  specifically,  the  SMSE  framework  enables  a  systematic 
approach  to  waveform  design  from  operations  research  -  a  field  of  study  dedicated  to  the 
various  forms  of  optimization  [18]. 

Coexistent  interference  occurs  when  two  or  more  communication  systems  operate 
without  orthogonality  in  frequency,  time,  space,  polarization,  and/or  coding.  With  limited 
available  bandwidth,  it  is  often  necessary  for  signals  to  spectrally  coincide  while  inducing 
“manageable”  levels  of  mutual  interference.  When  user  requirements  dictate  fundamentally 
different  waveform  modulations,  the  system  design  procedure  often  involves  trial  and  error 
to  find  waveforms  which  can  coexist  [25,  31].  As  commonly  employed  in  the  operations 
research  field,  the  Genetic  Algorithm  (GA)  and  Response  Surface  Methodology  (RSM) 
techniques,  as  described  in  Section  2.2.1  and  Section  2.2.2,  represent  two  approaches  that 
can  be  considered  for  waveform  design  optimization  to  provide  a  more  structured,  optimal 
means  for  determining  waveform  design  parameters  [3,4]. 

Each  of  these  optimization  techniques  are  employed  here  to  demonstrate  SMSE 
waveform  design  in  a  coexistent  scenario  containing  an  SMSE  signal  and  a  Direct  Se¬ 
quence  Spread  Spectrum  (DSSS)  system  operating  over  an  Additive  White  Gaussian  Noise 
(AWGN)  channel.  The  DSSS  system  employs  a  non-coherent  Delay-Eock  Eoop  (DEE) 
as  described  in  Section  2.1.3  for  code  tracking  prior  to  data  demodulation.  It  is  important 
to  note  that  the  particular  DEE  implementation  used  here,  and  the  metric  introduced  in 
Section  3.2.3  to  characterize  various  tracking  conditions  (perfect  and  imperfect),  are  suffi¬ 
ciently  general  such  that  the  optimization  demonstrations  herein  are  broadly  applicable  to 
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other  non-communication  applications  employing  DLL  tracking,  e.g.,  precision  navigation, 
timing,  geolocation,  etc. 

The  coexistent  SMSE-DSSS  waveform  design  process  is  demonstrated  under  two 
conditions,  including:  1)  perfect  DLL  code  tracking  as  described  in  Section  3.1  and  2)  im¬ 
perfect  DLL  code  tracking  as  described  in  Section  3.2.  Under  perfect  DLL  code  tracking 
conditions,  the  GA  and  RSM  optimization  processes  are  considered  independently  in  Sec¬ 
tion  3.1.3  and  Section  3.1.4,  respectively.  In  both  cases,  the  optimization  goal  is  to  deter¬ 
mine  SMSE  parameter  values  for  the  number  of  lEET  points  Nf  and  subcarrier  spacing  A/ 
such  that  DSSS  bit  error  rate  Pb  is  optimized  (minimized  or  maximized).  Under  imperfect 
DEE  code  tracking  conditions  in  Section  3.2,  the  GA  and  RSM  techniques  are  sequentially 
combined  into  a  hybrid  optimization  process  that  includes:  1)  the  GA  process  being  ap¬ 
plied  in  Section  3.2.4  to  generate  a  “coarse”  solution  for  initial  RSM  processing,  and  2)  the 
RSM  process  providing  the  final  optimized  solution  in  Section  3.2.5.  As  in  the  perfect  code 
tracking  case,  the  end  goal  is  to  determine  the  {N A/)  pair  that  optimizes  Pb-  However, 
Pb  optimization  (minimization  or  maximization)  is  actually  accomplished  through  a  Corre¬ 
lation  Degradation  metric  Coeg  as  introduced  in  Section  3.2.3.  The  statistical  behavior  of 
Coeg  is  used  to  capture  and  characterize  overall  DEE  code  tracking  performance. 

3.1  Perfect  DLL  Code  Tracking 

3.1.1  Coexistent  SMSE-DSSS  Scenario.  The  following  signal  conditions  were 
used  for  coexistent  SMSE-DSSS  scenario  under  perfect  DEE  code  tracking  conditions. 

The  SMSE  signal  was  generated  according  to  the  framework  described  in  Section  2.1.1. 
Two  SMSE  factors  (design  parameters)  were  varied  for  the  experiments,  including:  1)  the 
total  number  of  lEET  points  Nf  and  2)  the  subcarrier  frequency  separation  A/.  The  re¬ 
maining  SMSE  design  parameters  in  (2.4)  were  fixed  such  that  conventional  OEDM  was 
implemented  [34],  i.e.,  Sk  =  dk.  The  complex  baseband  OEDM  symbols  were  generated 
using  independent  BPSK  data  modulation  on  all  subcarriers  and  carrier  modulated  to  fc  for 
coexistent  demonstrations.  The  carrier  modulated  SMSE  waveform  occupies  a  total  band¬ 
width  of  WsMSE  =  2  X  Nf  X  A/  and  has  a  duration  of  Tqfdm  =  ^/Rofdm  =  1/A/. 
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The  SMSE  data  bits  were  randomly  generated  with  equal  probability.  Therefore, 
every  transmitted  SMSE  symbol  was  random  and  there  were  no  speeial  bit  sequenees  used 
to  simulate  frames,  paekets,  ete.  Einally,  there  was  no  eyelie  prefix  used  for  the  waveform 
design  demonstrations  and  analysis. 

The  eoexistent  DSSS  system  used  BPSK  for  both  data  and  spreading  modulations. 
Eor  demonstration  purposes,  the  spreading  eode  was  an  =  32-bit  Hadamard  sequenee 
with  exaetly  one  eode  period  (32  ehip  intervals)  oeeurring  per  data  symbol.  The  DSSS 
symbol  duration  was  Tsym  =  1  sec,  and  the  symbol  rate  was  Rsym  =  ^/Tsym  =  I  Hz. 
Consequently,  the  ehip  rate  was  Rchip  =  NcX  Rsym  =  32  Hz.  The  DSSS  reeeiver  was  per¬ 
fectly  synchronized  to  the  transmitted  DSSS  signal  in  terms  of  earner  traeking  (frequeney 
and  phase),  spreading  eode  traeking,  and  eommunieation  symbol  traeking.  An  ideal  RE  fil¬ 
ter  was  used  and  eommunieation  symbols  were  estimated  using  a  single  ehannel  eorrelation 
reeeiver  under  Maximum  Eikelihood  (ME)  eonditions. 

3.1.2  Optimization  Metric.  The  optimization  metrie  under  perfeet  DEE  eode 
traeking  traeking  eonditions  was  end-to-end  DSSS  system  Pj,,  as  determined  by  Monte 
Carlo  simulation  of  a  physies-based  analytie  model.  The  model  assumed  that  the  eoexistent 
SMSE  and  DSSS  signals  were  speetrally  eoineident  (same  eenter  frequeney)  and  were 
operating  over  an  AWGN  ehannel.  Therefore,  the  resultant  DSSS  bit  errors  are  due  to  a 
eombination  of  ehannel  noise  and  the  eoexistent  SMSE  signal.  However,  the  ehannel  noise 
power  was  fixed  during  both  the  minimization  and  maximization  demonstrations.  Thus,  the 
DSSS  Pb  eurves  in  Seetion  3.1  eorrespond  to  bit  error  ehange  as  a  funetion  of  interfering 
signal  power  and  optimized  input  SMSE  parameters.  The  ehanges  are  not  due  to  differing 
noise  power. 

3.1.3  Genetic  Algorithm  ( GA).  Eaeh  eombination  of  the  two  optimization  input 
parameters  (Nf,  A/)  were  represented  using  eight  binary  digits  in  a  gene.  Eurthermore, 
the  number  of  SMSE  subearriers  was  eonstrained  to  be  an  integer  power  of  two  with  Nf  G 
[1, 128].  Similarly,  the  SMSE  subearrier  spaeing  was  assigned  an  integer  value  satisfying 
A/  G  [1,  33].  The  feasibility  region  for  optimization  ineluded  all  possible  eombinations 
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of  Nf  and  A/  within  these  ranges.  For  the  special  case  of  Nf  =  1,  the  resultant  SMSE 
waveform  is  equivalent  to  a  single  BPSK  modulated  subcarrier  and  has  an  RF  bandwidth 
of  WsMSE  =  ‘^/Tqfdm  =  2/A/,  where  A/  is  used  here  and  in  presenting  subsequent 
Nf  =  1  results  to  more  appropriately  refer  to  waveform  bandwidth  and  not  subcarrier 
spacing.  Finally,  the  objective  function  was  end-to-end  DSSS  (P^). 

Consistent  with  the  process  described  in  Section  2.2.1,  the  GA  process  proceeded  as 
follows: 

1 .  The  initial  GA  population  consisted  of  ten  randomly  generated  genes. 

2.  Using  Monte  Carlo  simulation  of  a  physics-based  analytic  model,  the  DSSS  P^  was 
calculated  for  each  case  to  judge  each  genes’  fitness. 

3.  For  mating,  four  of  the  five  most  fit  genes  and  two  of  the  least  fit  genes  were  chosen 
as  parents.  The  resulting  six  parents  were  randomly  assigned  for  mating,  creating 
three  pairs  of  parents. 

4.  To  mate,  the  binary  digits  of  the  parents’  genes  were  compared.  When  the  binary  val¬ 
ues  of  parent  genes  matched,  the  same  values  were  passed  to  the  offspring.  Where 
the  values  differed,  the  offspring  values  were  randomly  chosen  with  equal  probabil¬ 
ity.  Each  pair  of  parents  created  two  offspring. 

5.  After  mating,  each  binary  digit  in  the  offspring’s  gene  was  subjected  to  a  10%  chance 
of  mutation,  or  complementing  the  bit  value. 

6.  The  next  population  of  ten  genes  included  the  six  children  and  the  four  most  fit 
parents.  The  GA  process  then  proceeded  as  it  did  with  the  initial/previous  population 
of  ten  genes.  For  design  demonstrations  in  this  dissertation,  the  GA  process  was 
repeated  for  100  generations. 

The  GA  optimization  process  was  used  to  both  minimize  and  maximize  DSSS  Pi,. 
Minimizing  Ph  creates  what  could  be  called  “peaceful”  coexistence  while  maximizing 
creates  worst-case  coexistence  conditions.  Results  for  GA  optimization  under  perfect  code 
tracking  conditions  are  provided  in  Section  4.1.2. 
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3.1.4  Response  Surface  Methodology  (RSM).  The  RSM  process  was  also  used  to 
optimize  SMSE  parameters  (Nj,  A/)  such  that  coexistent  DSSS  performance  was  both 
maximized  and  minimized.  The  RSM  experiments  assumed  a  second-order  model  as  de¬ 
scribed  in  (2.17).  As  a  result,  the  experimental  design  consisted  of  a  two-factor,  three-level, 
full-factorial  design  with  four  additional  center  runs.  The  corresponding  system  model 
from  (2.8)  for  this  demonstration  is  expressed  as 


1-1-1  111 
1-1  0  0  10 

1-1  1-111 
1  0-1  0  0  1 

1  0  0  0  0  0 

1  0  1  0  0  1 


Y  =  1  1-1-111 

1  1  0  0  10 


(3.1) 


1  0  0  0  0  0 
1  0  0  0  0  0 
1  0  0  0  0  0 
1  0  0  0  0  0 


The  SMSE  parameters  were  optimized  using  the  steepest  ascent/descent  process  described 
in  (2.21)  until  the  response  surface  fit  the  second-order  model,  as  determined  by  the  ANOVA. 
Once  the  second-order  model  detected  curvature,  (2.19)  was  used  to  solve  for  the  SMSE 
parameters  in  terms  of  their  coded  variables.  Results  for  the  RSM  optimization  process 
under  perfect  code  tracking  conditions  are  provided  in  Section  4.1.3. 

3.2  Imperfect  DLL  Code  Tracking 

3.2.1  Coexistence  Scenario.  The  following  signal  conditions  were  used  for  co¬ 
existent  SMSE-DSSS  scenario  under  imperfect  DEE  code  tracking  conditions. 
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The  SMSE  coexistent  signal  remained  unchanged  from  that  used  in  the  perfect  DLL 
code  tracking  demonstration  in  Section  3.1.1.  Once  gain,  the  DSSS  system  used  BPSK 
for  both  data  and  spreading  modulations.  However,  the  spreading  code  used  for  demon¬ 
stration  was  an  Nc  =  1023-chip  Gold  code  sequence  with  exactly  one  code  period  (1023 
chip  intervals)  occurring  per  data  symbol.  The  symbol  rate  was  set  to  Rsym  =  1  kHz, 
and  therefore  the  chip  rate  was  Rchip  =  1.023  MHz.  The  DSSS  receiver  was  perfectly 
synchronized  to  the  transmitted  DSSS  signal  in  terms  of  carrier  tracking  (frequency  and 
phase)  and  communication  symbol  tracking. 

The  DSSS  receiver  used  the  non-coherent  DLL  described  in  Section  2.1.3  for  code 
tracking.  A  non-coherent  DLL  was  chosen  for  demonstration  give  it  is  suitable  for  gen¬ 
eral  purpose  DSSS  receivers  without  requiring  precise  carrier  tracking  [25].  The  DLL  RL 
filter  was  an  8*^-order  Butterworth  filter  having  a  bandwidth  equal  to  twice  the  chip  rate, 
Wrf  =  2Rchip  =  2.046  MHz.  The  bandpass  filter  following  the  despreading  mixer  in 
the  early/late  DLL  branches  was  an  8*^-order  Butterworth  filter  having  a  bandwidth  of 
Wbpf  =  5  kHz.  The  low-pass  filter  in  the  DLL  energy  detectors  used  a  4*^-order  Butter¬ 
worth  filter  with  a  bandwidth  of  Wlpf  =  2.5  MHz.  The  loop  filter  was  a  first-order  filter 
with  Fioop{s)  =  1.  The  mapping  between  the  discriminator  output  and  the  NCO  was  linear, 
such  that  a  maximum  response  from  the  discriminator  resulted  in  a  NCO  code  generation 
rate  of  2Rchip-  Linally,  the  early  (advanced)  and  late  (delayed)  codes  were  separated  by 
Ae-i  =  1  chip. 
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3.2.2  Hybrid  Optimization.  The  perfect  DLL  code  tracking  demonstration  in 
Section  3.1  only  accounted  for  a  portion  of  the  errors  associated  with  coexistent  interfer¬ 
ence.  Under  more  realistic  conditions,  imperfect  DLL  code  tracking  due  to  interfering 
signals  will  further  degrade  DSSS  Pb  performance  [15, 16].  The  imperfect  code  tracking 
demonstration  is  designed  to  help  isolate  code  tracking  error  effects  resulting  from  a  coex¬ 
istent  SMSE  waveform  being  received  by  the  DSSS  system. 

Unlike  the  perfect  DLL  code  tracking  demonstration  in  Section  3.1,  the  GA  and 
RSM  optimization  processes  are  not  performed  independently  under  imperfect  tracking 
conditions.  Rather,  they  are  used  to  complement  each  other  using  a  hybrid  optimization 
approach.  In  this  process,  the  GA  process  is  first  used  to  determine  an  initial  “coarse”  solu¬ 
tion  that  is  then  passed  to  the  RSM  process  which  determines  the  final  optimized  solution. 
This  approach  is  advantageous  for  several  reasons.  First,  given  that  GA  is  naturally  a  dis¬ 
crete  process,  it  performs  best  at  describing  and  optimizing  SMSE  design  parameters  that 
are  discrete  as  well,  e.g.,  the  number  of  IFFT  points  Nf.  Second,  GA  solutions  are  less 
accurate  when  the  gene  mapping  is  for  continuous  variables,  such  as  subcarrier  spacing 
A/.  In  this  respect,  the  RSM  process  is  most  advantageous  given  its  solutions  are  based 
on  a  modeled  surface  response  and  the  resultant  optimized  solution  is  not  necessarily  part 
of  the  input  test  matrix.  However,  one  limitation  of  the  RSM  process  is  that  it  requires  a 
good  starting  point  or  the  search  process  may  become  too  prolonged.  Therefore,  the  final 
hybrid  approach  for  optimizing  (Nf,Af)  selection  exploits  the  strength  of  each  process 
and  consists  of  1)  using  GA  first  to  determine  the  most  appropriate  Nf  value,  followed  by 
2)  the  RSM  process  to  find  the  optimized  A/  value  associated  with  the  GA  Nf  value. 
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3.2.3  Optimization  Metric.  To  isolate  the  impact  of  a  coexistent  SMSE  signal 
on  DSSS  code  tracking  performance,  an  optimization  metric  is  introduced.  In  this  case, 
optimization  is  based  on  DLL  code  tracking  performance  versus  end-to-end  DSSS  Pi,.  The 
DSSS  Correlation  Degradation  metric  is  defined  here  as 


(3.2) 


where  Cp{t)  is  the  DLL  prompt  code  estimate  shown  in  Ligure  2.10,  Ci{t)  is  incident  re¬ 
ceived  code  and  integration  is  carried  out  over  one  full  code  period  (N^Tc  =  Tgym  = 
1  msec  for  this  demonstration).  The  minimum  value  of  Coeg  =  0  indicates  perfect  code 
tracking,  a  value  of  0  <  Coeg  <  1  indicates  imperfect  code  tracking,  and  a  maximum 
value  of  Coeg  =  1  indicates  a  DLL  break-lock  condition.  Representative  histograms  for 
Coeg  in  (3.2)  under  imperfect  manageable  and  imperfect  severely  degraded  DLL  tracking 
conditions  are  shown  in  Lig.  3.1  and  Lig.  3.2,  respectively. 

3.2.4  Genetic  Algorithm  ( GA).  Each  combination  of  the  two  optimization  input 
parameters  (Nf,  A/)  were  represented  using  nine  binary  digits  in  a  gene.  In  addition, 
the  parameters  were  constrained  to  Nf  G  {16,  32,  64, 128}  and  A/  e  [6, 133]  fl  X.  The 
feasibility  region  was  defined  such  that  the  total  SMSE  signal  bandwidth  was  less  than  the 
DLL  RL  filter  bandwidth  (NfAf  <  2.046  MHz).  The  optimization  objective  function  was 
the  the  correlation  metric  Coeg  given  in  (3.2). 

Consistent  with  the  process  described  in  Section  2.2.1,  the  GA  process  proceeded  as 
follows: 
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Correlation  Degredation,  Coeg 


Figure  3.1:  Representative  histogram  of  eorrelation  degradation  metrie 
Coeg  in  (3.2)  for  imperfeet  manageable  DLL  traeking  eonditions. 


Correlation  Degredation,  Coeg 


Figure  3.2:  Representative  histogram  of  eorrelation  degradation  metrie 
Coeg  in  (3.2)  for  imperfeet  severely  degraded  traeking  eonditions. 
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1.  The  initial  population  consisted  of  ten  randomly  generated  genes. 

2.  Using  a  physics-based  analytic  model,  Coeg  was  calculated  for  each  case  to  judge 
the  genes’  fitness. 

3.  For  mating,  four  of  the  five  most  fit  genes  and  two  of  the  least  fit  genes  were  chosen 
as  parents.  The  resulting  six  parents  were  randomly  assigned  for  mating,  creating 
three  pairs  of  parents. 

4.  To  mate,  the  binary  digits  of  the  parents’  genes  were  compared.  When  the  binary  val¬ 
ues  of  parent  genes  matched,  the  same  values  were  passed  to  the  offspring.  Where 
the  values  differed,  the  offspring  values  were  randomly  chosen  with  equal  probabil¬ 
ity.  Each  pair  of  parents  created  two  offspring. 

5.  After  mating,  each  binary  digit  in  the  offspring’s  gene  was  subjected  to  a  10%  chance 
of  mutation,  or  complementing  the  bit  value. 

6.  The  next  population  of  ten  genes  included  the  six  children  and  the  four  most  fit 
parents.  The  GA  process  then  proceeded  as  it  did  with  the  initial/previous  population 
of  ten  genes.  For  design  demonstrations  in  this  dissertation,  the  GA  process  was 
repeated  for  1000  generations. 

The  GA  process  was  used  to  both  minimize  and  maximize  the  DSSS  Coeg-  Minimiz¬ 
ing  Coeg  corresponds  to  “peaceful”  coexistence  while  maximizing  Coeg  creates  a  worst- 
case  coexistence  scenario.  The  GA  output  included  ten  SMSE  {N f,  A/)  parameter  com¬ 
binations  for  both  the  minimization  and  maximization  cases.  These  final  populations  were 
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compared  using  multi-comparison  tests  from  Section  2. 2. 2. 8.  Optimization  results  for  this 


process  are  found  in  Section  3.2.4. 


3.2.5  Response  Surface  Methodology  (RSM).  Given  that  the  GA  process  pro¬ 
vided  a  “coarse”  solution  that  maximized  and  minimized  Coeg,  the  RSM  was  next  used  to 
obtain  final  optimization  results.  Since  Nj  is  a  discrete  parameter  and  every  value  is  tested 
in  the  GA  process,  the  GA-optimized  values  for  Nf  were  accepted.  The  RSM  process  was 
then  applied  to  find  corresponding  optimized  values  for  A/. 


As  in  the  perfect  code  tracking  case,  the  RSM  experiments  used  the  ANOVA  with  a 
second-order  model  to  detect  curvature  in  the  response  surface.  However,  there  was  only 
one  input  variable  in  the  model,  xi.  Consequently,  the  experimental  design  consisted  of  a 
one-factor,  five-level,  full-factorial  design.  The  matrix  form  of  the  system  model  from  (2.8) 
is  then  expressed  as 


Y  = 


1  -2  4 
1  -1  1 
1  0  0 
1  1  1 
1  2  4 


/^o 

Pi 

Pii 


(3.3) 


The  SMSE  parameters  were  optimized  using  the  steepest  ascent/descent  process  de¬ 
scribed  in  (2.21)  until  the  response  surface  fit  the  second-order  model,  as  determined  by  the 
ANOVA.  Once  the  second-order  model  detected  curvature,  (2.19)  was  used  to  solve  for  the 
SMSE  parameters  in  terms  of  their  coded  variables.  Einal  results  for  the  RSM  optimization 
process  under  imperfect  code  tracking  conditions  are  provided  in  Section  4.2.3. 
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IV.  Results  and  Analysis 

The  SMSE  waveform  design  proeedure  is  demonstrated  in  this  ehapter  under  two  code 
tracking  conditions,  including  perfect  code  tracking  and  imperfect  code  tracking.  In  both 
cases,  the  SMSE  waveform  is  introduced  into  the  DSSS  system  and  resultant  DSSS  bit  error 
rate  (P^)  performance  characterized  as  the  SMSE  parameters  are  varied.  Under  perfect  code 
tracking  conditions  in  Section  4.1,  the  DSSS  receiver  is  assumed  to  maintain  perfect  code 
tracking  such  that  the  DEE  produces  an  ideal  prompt  code  Cp{t)  for  despreading.  Optimal 
SMSE  parameter  selection  is  addressed  using  independent  Genetic  Algorithm  (GA)  and 
Response  Surface  Methodology  (RSM)  approaches.  Under  imperfect  code  tracking  condi¬ 
tions  in  Section  4.2,  DEE  code  tracking  is  incorporated  to  assess  performance  under  more 
realistic  channel  conditions  using  a  less  than  ideal  Cp{t)  estimate.  Optimal  SMSE  parame¬ 
ter  selection  under  these  conditions  is  addressed  using  a  hybrid  GA-RSM  technique.  The 
reader  should  exercise  caution  when  comparing  optimization  results  across  various  sec¬ 
tions  in  this  chapter.  This  is  especially  important  when  considering  perfect  and  imperfect 
DSSS  code  tracking  results,  given  they  were  generated  using  dissimilar  SMSE  and  DSSS 
parameter  values  as  well  as  different  objective  functions  and  feasibility  region  constraints 
during  optimization. 

4.1  Perfect  Code  Tracking 

4.1.1  Demonstration  Procedure.  The  SMSE  waveform  design  procedure  is  first 
demonstrated  in  a  coexistent  environment  under  perfect  DSSS  code  tracking  conditions. 
In  this  case,  the  estimated  DEE  prompt  code  Cp{t)  is  considered  to  be  ideal  such  that  the 
DSSS  receiver  despreading  code  perfectly  matches  the  transmitted  spreading  code.  Opti¬ 
mal  SMSE  parameter  selection  is  addressed  using  independent  GA  and  RSM  approaches. 
The  goal  is  to  find  SMSE  parameter  values  that  optimize  DSSS  receiver  performance  in 
terms  of  end-to-end  bit  error  rate  (Pfe). 

4. 1.1.1  Coexistent  SMSE  Signal.  The  coexisting  SMSE  signal  was  gen¬ 
erated  using  the  framework  described  in  Section  2.1.1  all  but  two  of  the  parameters  fixed 
to  implement  conventional  OEDM  [34].  Performance  of  an  SMSE  OEDM  implementa- 
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tion  is  dictated  by  two  design  parameters,  ineluding  the  total  number  of  IFFT  points  Nf, 
whieh  dietates  the  number  of  individual  SMSE  subearriers,  and  the  subcarrier  frequeney 
spaeing  A/.  Together,  the  value  of  these  two  parameters  determines  the  overall  waveform 
bandwidth. 


4. 1.1. 2  Overall  DSSS  System.  The  DSSS  system  eonsidered  here  uses 
BPSK  for  both  data  and  spreading  modulations.  The  spreading  eode  was  a  iVc  =  32  bit 
Hadamard  sequenee  with  exaetly  one  eode  period  oeeurring  per  data  symbol.  The  DSSS 
Pb  performanee  is  evaluated  over  an  Additive  White  Gaussian  Noise  (AWGN)  ehannel. 
The  DSSS  reeeiver  employs  an  ideal  RF  filter  prior  to  the  despreading  mixer.  Following 
the  despreading  mixer,  the  eommunieation  symbols  are  estimated  using  a  single  ehannel 
eorrelation  proeess  under  maximum  likelihood  eonditions. 

4.1.2  Genetic  Algorithm  ( GA ).  The  analytie  SMSE  framework  enables  paramet- 
rie  optimization  of  OEDM  waveforms  using  a  GA  proeess.  As  noted  previously,  the  two 
SMSE  optimization  parameters  ineluded,  1)  the  number  of  IFFT  points  in  the  SMSE  wave¬ 
form  generation  (Nf),  and  2)  the  subearrier  frequeney  separation  (A/).  Eaeh  GA  gene 
used  eight  binary  digits  to  represent  the  possible  SMSE  parameter  values.  Eurthermore, 
the  value  of  Nf  was  eonstrained  to  be  an  integer  power  of  two  in  the  range  [1,128].  Eor 
representation  in  the  GA  gene.  A/  was  assigned  an  integer  value  in  the  range  [1,33].  The 
objeetive  funetion  to  be  optimized  was  DSSS  bit  error  rate  Pb-  Eor  demonstration  purposes, 
the  GA  optimization  proeess  was  earned  out  for  two  eases:  (A)  Minimizing  DSSS  Pb  whieh 
represents  best-ease  SMSE-DSSS  eoexistenee  performanee,  and  (B)  Maximizing  DSSS  Pb 
whieh  represents  worst-ease  SMSE-DSSS  eoexistenee  performanee. 

Eor  both  oases,  the  initial  GA  population  oonsisted  of  ten  randomly  generated  genes. 
Using  Monte  Carlo  simulation  of  a  physios-based  analytie  model,  the  DSSS  Pb  was  oal- 
oulated  for  eaeh  ease  to  judge  the  genes’  fitness,  with  the  DSSS  system  Pb  value  being 
the  fitness  statistie.  Tabulated  GA  results  for  the  two  optimization  eases  are  presented  in 
Table  4.1  and  Table  4.2  [3].  Eaeh  table  inoludes  the  initial  random  population  values  and 
the  final  optimized  values.  Table  4.1  shows  the  GA  optimized  parameter  values  that  min- 
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imize  DSSS  Pf,  (best-case  SMSE-DSSS  coexistence)  and  Table  4.2  shows  GA  optimized 
parameter  values  that  maximize  DSSS  Pb  (worst-case  SMSE-DSSS  coexistence). 

Corresponding  bit  error  curves  for  the  data  presented  in  Table  4.1  and  Table  4.2  are 
shown  in  Eig.  4.1  and  Eig.  4.2,  respectively  [3].  In  both  cases,  the  signal  power  (S)  and 
the  noise  power  (N)  in  the  DSSS  system  remained  fixed  while  the  interfering  SMSE  power 
(/)  was  varied  to  achieve  the  indicated  Signal-to-Interference-plus-Noise  Ratio  SINK  = 

S/il  +  N). 

GA  optimized  results  for  Minimum  DSSS  Pb  are  shown  in  Eig.  4.1.  The  upper 
curve  (unfilled  boxes)  is  provided  for  comparison  and  shows  the  highest  resultant  Pb  per¬ 
formance  obtained  from  the  initial  non-optimized  parameter  population.  The  lower  curve 
(filled  boxes)  is  the  resultant  Pb  performance  using  the  final  GA-optimized  parameters  from 
Table  4.1  (Nj  =  64,  A/  =  11  Hz,  Wsmse  =  1-41  KHz)  and  represents  best-case  coex¬ 
istence.  In  this  case,  the  Wsmse  bandwidth  greatly  exceeds  Wrf  =  2  x  Rchip  =  Q4:  Hz 
and  minimal  received  SMSE  power  actually  enters  the  DSSS  detector.  Consequently,  the 
Pb  results  approach  the  noise  limited  case,  i.e.,  the  resultant  Pb  is  the  same  as  if  there  were 
no  SMSE  signal  present.  This  result  is  a  direct  consequence  of  the  feasibility  region  being 
loosely  constrained  and  allowing  solutions  where  Wsmse  >  Wrf-  As  such,  there  are  ad¬ 
ditional  GA  solutions  in  Table  4.1  that  were  analyzed  and  produced  similar  noise  limited 
results,  e.g.,  the  (Nj  =  128,  A/  =  11  Hz,  Wsmse  =  2.82  KHz)  solution  uses  the  same 
subcarrier  spacing  with  more  carriers  and  an  even  wider  bandwidth. 

GA  optimized  results  for  Maximum  DSSS  Pb  degradation  are  shown  in  Eig.  4.2.  The 
lower  curve  (unfilled  boxes)  is  provided  for  comparison  and  shows  the  lowest  resultant  Pb 
performance  obtained  from  the  initial  non-optimized  parameter  population  (noise  limited 
performance).  The  upper  curve  (filled  boxes)  is  the  resultant  Pb  performance  using  the  final 
GA-optimized  parameters  from  Table  4.2  (Nf  =  1,  /  =  16  Hz,  Wsmse  =  32  Hz).  This  is 
the  special  case  of  Nf  =  1  (single  BPSK  modulated  subcarrier  centered  at  fc)  and  produces 
worst-case  coexistence.  Given  the  resultant  bandwidth  of  Wsmse  =  ^2  Hz,  all  SMSE 
power  is  within  Wrf  =  64  Hz  and  contributes  to  degraded  performance.  In  addition  to 
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Table  4.1:  Initial  and  optimized  SMSE  parametric  value  populations  us¬ 
ing  the  GA  process  to  MINIMIZE  DSSS  probability  of  bit  error  {P^)  -  Best- 
Case  SMSE-DSSS  Coexistence  Performance  [3]. 


Nf 

A/  (kHz) 

Initial 

Optimized 

Initial 

Optimized 

128 

64 

20 

10 

2 

128 

28 

11 

16 

128 

28 

4 

8 

64 

29 

22 

128 

8 

5 

1 

64 

64 

4 

8 

8 

128 

21 

11 

1 

128 

5 

11 

64 

64 

19 

23 

8 

64 

18 

11 

Table  4.2:  Initial  and  optimized  SMSE  parametric  value  populations  us¬ 
ing  the  GA  process  to  MAXIMIZE  DSSS  probability  of  bit  error  (Pfe)  - 
Worst-Case  SMSE-DSSS  Coexistence  Performance  [3]. 


Nf 

A/  (kHz) 

Initial 

Optimized 

Initial 

Optimized 

128 

1 

20 

16 

2 

4 

28 

16 

168 

32 

28 

13 

8 

1 

29 

16 

128 

1 

5 

32 

64 

1 

4 

16 

8 

1 

21 

16 

1 

1 

5 

16 

64 

1 

19 

16 

8 

1 

18 

16 
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Signal-to-Interference  Plus  Noise  Ratio,  SINK  (dB) 


Figure  4.1:  SINK  vs  DSSS  Pb  performance  for  initial  (unfilled  boxes)  and 
GA  optimized  (filled  boxes)  SMSE  parameters  in  Table  4.1.  GA  optimized 
fox  Minimum  Pb  -  Best-Case  SMSE-DSSS  Coexistence  Performance  [3]. 
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Signal-to-Interference  Plus  Noise  Ratio,  SINR  (dB) 


Eigure  4.2:  SINR  vs  DSSS  Pb  performance  for  initial  (unfilled  boxes)  and 
GA  optimized  (filled  boxes)  SMSE  parameters  in  Table  4.2.  GA  optimized 
for  Maximum  Pb  -  Worst-Case  SMSE-DSSS  Coexistence  Performance  [3]. 
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power,  the  Pb  degradation  is  also  attributable  to  the  spectral  sine  (/)  =  sin  (/)  //response 
of  the  SMSE  waveform  and  the  spectral  structure  of  the  DSSS  spreading  signal.  More 
specifically,  the  discrete  spectral  lines  in  the  Nc  =  32-chip  DSSS  spreading  code  are  spaced 
fchip  =  Rchip/Nc  =  I  Hz  apart  [25]  and  there  are  WRp/fchip  =  64  total  lines  in  Wrf  = 
2  X  Rchip  =  64  Hz,  or  32  total  lines  in  Wsmse  =  Wrf/2  =  32  Hz.  The  impact  of  this 
is  illustrated  by  considering  the  power  spectrum  at  the  despreading  mixer  output  (DSSS 
detector  input)  in  the  DSSS  receiver,  which  is  the  convolution  of  the  spreading  code  spectral 
lines  and  the  SMSE  subcarriers.  Eor  the  shift-multiply-integrate  operations  of  the  spectral 
convolution  process  that  occur  near  /„  there  are  32  contributing  products  in  the  correlation 
result.  As  a  result,  the  amount  of  power  in  the  “despread”  spectral  response  that  falls  within 
the  DSSS  detector  bandwidth  (lEoet  =  2  x  Rsym  =  2  Hz)  is  maximum  and  the  interfering 
SMSE  signal  has  maximum  impact  on  DSSS  Pb  performance. 

4.1.3  Response  Surface  Methodology  (RSM).  The  RSM  was  the  second  approach 
considered  to  optimize  SMSE  parameter  selection  under  perfect  DSSS  code  tracking  con¬ 
ditions.  In  this  case,  the  experiments  assumed  the  Pb  response  surface  fit  a  second-order 
model  given  by 


Y  —  /3o  +  /3lXl  -f  [32^2  +  /3i2XiX2  -f  /3llXii  -|-  (322^22  (4.1) 

To  properly  characterize  this  model,  the  experimental  design  considered  a  two-factor,  three- 
level,  full-factorial  design  with  four  additional  center  runs.  The  matrix  form  of  the  experi¬ 
mental  design  is  given  by 
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Y 


1-1-1  111 
1-1  0  0  10 

1-1  1-111 
1  0-1  0  0  1 

1  0  0  0  0  0 

1  0  1  0  0  1 

1  1-1-111 
1  1  0  0  10 

1  1  1  111 
0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 

0  0  0  0  0  0 


f^o 

A 

A 

A2 

Ai 

A2 


(4.2) 


The  mapping  from  SMSE  parameter  variables  (iVj,  A/)  to  eoded  variables  {xi,X2) 
for  the  initial  and  final  RSM  experiments  is  shown  in  Table  4.3  and  Table  4.4,  respee- 
tively  [4].  The  final  results  were  obtained  using  the  method  of  steepest  deseent. 


Table  4.3:  Mapping  from  SMSE  parameters  (Aj-,  A/)  to  coded  variables 
(xi,  X2)  for  the  initial  RSM  experiment  [4]. 


Eactor  Eevel 

Nf 

Xi 

A/ 

X2 

Eow 

2^ 

-1 

17 

-1 

Medium 

2^ 

0 

16 

0 

High 

2^ 

1 

15 

1 

Table  4.4:  Mapping  from  SMSE  parameters  (A/,  A/)  to  coded  variables 
(xi,  X2)  for  the  final  RSM  experiment  after  applying  the  method  of  steepest 
descent [4]. 


Eactor  Eevel 

Nf 

Xi 

A/ 

X2 

Eow 

2^ 

-1 

16.37 

-1 

Medium 

22 

0 

15.87 

0 

High 

2^ 

1 

15.36 

1 
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Table  4.5:  ANOVA  table  for  initial  and  final  RSM  trials  [4] 


Souree  of 
Variation 

Degrees  of 
Freedom 

Initial 

p-value 

Final 

p-value 

Model 

5 

3.8  X  10-6 

2.1  X  10-^ 

Xi 

1 

1.6  X  10-® 

1.0  X  10-6 

X2 

1 

1.0  X  10“^ 

1.4  X  10-1 

X1X2 

1 

3.1  X  10“^ 

2.5  X  10-2 

Xi 

1 

6.6  X  10“^ 

1.5  X  10-6 

^2 

1 

3.8  X  10-1 

1.2  X  10-2 

Error 

7 

Table  4.5  shows  ANOVA  results  for  the  initial  and  final  RSM  experiments  [4],  with 
the  table  rows  representing  main  effeets  and  interaetions  aeeording  to  the  experimental 
model.  The  p-value  indieates  the  signifieanee  for  eaeh  faetor.  The  smaller  the  p-value, 
in  a  given  row,  the  more  likely  that  term  is  signifieant.  For  the  trials  illustrated  here,  a 
p-value  <  0.05  was  eonsidered  signifieant,  and  results  in  rejeetion  of  the  ANOVA  null 
hypothesis  of  equal  means. 

As  one  may  observe,  results  for  first  RSM  trial  indieate  that  the  response  surfaee  was 
not  quadratie  given  the  p- value  >  0.05.  After  moving  to  the  final  region,  both  of  the  xf 
and  X2  quadratie  terms  are  signifieant  and  the  stationary  point  Xg,  found  using  (2.19),  is 
determined  to  be  Nf  =  2  and  A/  =  15.87  Hz.  However,  additional  Eigenvalue  analysis 
of  the  resulting  regression  eoeffieients  revealed  that  the  stationary  point  Xg  is  neither  a 
minimum  nor  a  maximum,  but  rather  a  saddle  point.  If  a  global  minimum  or  maximum  is 
desired,  the  RSM  proeess  would  need  to  be  repeated  using  a  different  starting  point.  For 
purposes  of  this  researeh,  the  saddle  point  solution  is  suffieient  to  demonstrate  the  praetieal 
utility  of  the  RSM  proeess. 

To  illustrate  eonsisteney  between  the  RSM  saddle  point  solution  and  physieal  wave¬ 
form  level  modeling,  an  end-to-end  simulation  was  run  for  the  SMSE-DSSS  eoexistent 
seenario.  Simulation  results  are  shown  in  Fig.  4.3  for  the  Pb  eorresponding  to  the  RSM 
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Signal-to-Interference  Plus  Noise  Ratio,  SINK  (dB) 


Figure  4.3:  SINK  vs  DSSS  Pt,  performanee  for  RSM  optimized  SMSE 
parameters.  Unfilled  boxes  represent  Best-Case  SMSE-DSSS  Coexistence 
and  filled  boxes  represent  Worst-Case  SMSE-DSSS  Coexistence  [4]. 


stationary  point  (unfilled  circles)  [4].  Additional  best-case  (unfilled  boxes)  and  worst-case 
(filled  boxes)  results  are  also  shown  and  were  obtained  from  additional  RSM  searches. 

As  indicated  in  Eig.  4.3,  the  saddle  point  solution  in  Table  4.4  of  (Nf  =  4,  A/  = 
15.87  Hz,Wsmse  =  127  Hz)  has  approximately  one-half  of  its  power  within  Wrf  = 
2  X  Rchip  =  Q4:  Hz  and  is  indeed  non-optimal  (neither  maximum  nor  minimum  degrada¬ 
tion)  given  that  its  resultant  Pf,  performance  consistently  fall  between  the  two  extremes  of 
maximum  degradation  for  (Nf  =  1,  A/  =  8  Hz,  Wsmse  =  16  Hz)  and  minimum  noise 
limited  performance  for  (Nf  =  8,  A/  =  2  Hz,  Wsmse  =  32  Hz).  As  with  GA  maximiza¬ 
tion  results  in  Section  4.1.2,  worst-case  coexistence  performance  is  once  again  achieved 
for  the  special  case  of  Nf  =  1  (single  BPSK  modulated  subcarrier  centered  at  /c)  and  the 
physical  interpretation  as  to  why  this  occurs  is  as  explained  in  that  section. 

Results  in  Section  4.1.2  and  Section  4.1.3  suggest  that  the  GA  and  RSM  techniques 
are  applicable  for  rigorous  coexistence  analysis  of  conventional,  DSSS  and  OEDM-based 
SMSE  waveforms.  In  applying  both  techniques,  the  independently  optimized  results  con- 
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Figure  4.4:  Imperfect  code  tracking  scenario  for  SMSE-DSSS  coexis¬ 
tent  optimization.  SMSE  parameters  impact  DEE  performance  and  impact 
prompt  code  Cp{t)  estimation. 

sistently  migrated  toward  expected  system  performance  and  there  were  valuable  “lessons 
learned”  regarding  how  to  best  address  more  realistic  constraints  and  conditions.  As  pre¬ 
sented  in  the  next  section,  these  lessons  learned  drove  the  development  a  hybrid  GA-RSM 
optimization  process  to  address  imperfect  code  tracking  conditions. 

4.2  Imperfect  Code  Tracking 

The  SMSE  waveform  design  procedure  is  next  demonstrated  in  a  coexistent  envi¬ 
ronment  under  imperfect  code  tracking  conditions.  As  discussed  in  Section  2.1.3,  receiver 
code  tracking  is  critical  to  successful  demodulation.  The  effects  of  code  tracking  are  incor¬ 
porated  here  to  assess  SMSE-DSSS  coexistence  performance  under  more  realistic  channel 
conditions.  In  this  case,  a  less  than  ideal  prompt  code  estimate  Cp{t)  is  used  and  optimal 
SMSE  parameter  selection  is  addressed  using  a  hybrid  GA-RSM  technique.  The  hybrid 
technique  first  uses  the  GA  process  in  Section  4.2.2  to  find  a  “coarse”  optimization  solu¬ 
tion.  The  course  GA  solution  is  then  used  as  the  initial  starting  point  in  the  RSM  process 
of  Section  4.2.3  which  provides  the  final,  more  precise  optimized  solution. 

4. 2. 1  Demonstration  Procedure.  The  imperfect  code  tracking  scenario  for  SMSE- 
DSSS  optimization  is  depicted  in  Eig.  4.4  which  shows  the  coexistent  SMSE  and  DSSS 
transmitters.  As  shown,  both  signals  are  present  while  the  DEE  in  the  DSSS  receiver  tracks 
the  received  spreading  code. 
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Figure  4.5:  Coexistent  SMSE  signal  generation  arehiteeture. 

The  DLL  prompt  eode  output  Cp{t)  in  Lig.  4.4  represents  the  PN  spreading  eode 
estimate  that  is  used  to  despread  the  reeeived  signal  prior  to  demodulation.  In  this  see- 
nario,  the  reeeived  DSSS  waveform  is  eorrupted  by  the  SMSL  waveform  whieh  eauses  the 
prompt  eode  estimate  to  be  an  imperfeet  representation  of  the  true  spreading  eode.  As  a 
result,  the  ineoming  DSSS  signal  of  interest  is  not  perfeetly  despread  and  demodulated 
Pb  performanee  is  degraded.  Demonstration  details  with  respeet  to  the  eoexistent  SMSL 
signal,  the  overall  DSSS  system,  and  DLL  eode  traeking  are  presented  in  Seetion  4. 2. 1.1, 
Seetion  4. 2. 1.2  and  Seetion  4.2. 1.3,  respeetively. 

4.2. 1.1  Coexistent  SMSE  Signal.  The  SMSL  signal  was  generated  using 
the  arehiteeture  shown  in  Lig.  4.5.  The  sequenee  of  input  data  bits  dk  are  randomly  gen¬ 
erated  with  equal  probability  of  being  either  a  0  or  1 .  The  bits  are  then  modulated  using 
antipodal  BPSK  sueh  that  the  bit  is  mapped  to  dk  G  {  —  1, 1}  aeeording  to 

dk  =  (-1)'"=  . 

Using  a  serial-to-parallel  (S/P)  eonversion  proeess,  groups  of  Nf  modulated  bits 
passed  to  the  ILLT  operation  after  1)  element-by-element  weighting  by  eomplex  veetor 
Wk  and  2)  zero  padding.  The  first  weighting  eoeffieient  is  set  to  wi  =  0  and  the  remain¬ 
ing  Nf  —  1  eoeffieients  are  set  to  unity  sueh  that  Wk  =  [Oil  ■  ■  ■  1] .  This  is  eonsistent 
with  eommon  praetiee  when  implementing  OLDM  and  effeetively  ensures  that  the  resul- 
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tant  SMSE  waveform  does  not  include  a  DC  component.  In  the  simulation,  the  resultant 
weighted  vector  of  Nj  elements  is  then  zero  padded  with  a  total  of  fs/^f  —  Nf  zeros. 
This  ensures  that  the  resultant  time  domain  SMSE  symbol  following  the  lEET  process  is 
sampled  at  the  same  effective  rate  /*  as  other  signals  in  the  simulation.  In  the  simulation, 
this  operation  effectively  replaces  the  Digital-to-Analog  converter  (DAC)  in  a  real-world 
communication  system. 

Each  lEET  operation  creates  one  time  domain  SMSE  symbol  comprised  of  Nf  carri¬ 
ers  and  having  a  symbol  duration  of  Tsym  =  1/^/  and  sample  frequency  fs-  The  resultant 
SMSE  symbols  are  complex  baseband  signals.  The  received  interfering  SMSE  waveform 
was  taken  as  real  part  of  the  carrier  modulated  signal,  expressed  as 


ssMSE{t)  =  Re  , 

where  fc  is  the  carrier  frequency,  v{t)  is  the  complex  baseband  SMSE  signal,  and  ssMSsit) 
is  the  resultant  coexisting  SMSE  signal. 

The  optimized  SMSE  variables  used  for  the  imperfect  code  tracking  demonstration 
included:  1)  the  number  of  lEET  points  Nf  and  2)  the  subcarrier  frequency  separation  A/, 
which  were  allowed  to  take  on  values  of  Nf  G  {16,  32, 64, 128}  and  A/  G  [6, 132]  kHz. 

4.2. 1.2  Overall  DSSS  System.  The  transmitted  DSSS  signal  was  generated 
using  the  architecture  shown  in  Eig.  4.6.  The  input  data  modulated  waveform  d{t)  is  based 
on  BPSK  modulation  using  randomly  generated  bits  having  equal  probability  of  being  a 
0  or  1.  The  baseband  data  modulated  waveform  d{t)  is  then  carrier  modulated  to  fc  and 
spread  by  c{t)  prior  to  transmission.  Spreading  waveform  c{t)  is  generated  according  to 

c(f)  =  -  mTchip)  , 

m=l 

where  is  an  A}  length  binary  spreading  code  and  Tchip  is  the  chip  duration.  Eor  the 
imperfect  code  tracking  demonstration  being  considered  here:  1)  Cm  was  an  Nc  =  1023 
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Figure  4.6:  DSSS  signal  generation  arehiteeture. 

length  Gold  eode,  2)  the  eode  period  was  set  equal  to  the  symbol  duration  {Nc  x  Tchip  = 
Tsym),  3)  the  ehip  rate  was  Rchip  =  ^ /Tchip  =  1023  x  10^  ehips/see,  and  4)  the  eode  period 
start/stop  boundaries  were  synehronous  with  Tgym  transition  boundaries. 


4.2. 1.3  Delay-Lock  Loop  (DLL)  Code  Tracking.  The  DLL  arehiteeture 


used  for  the  imperfeet  eode  traeking  demonstration  is  shown  in  Fig.  4.7.  The  RF  bandpass 
filter  was  implemented  as  an  8*^-order  Butterworth  filter  having  a  —3  dB  bandwidth  of 

Wrf  =  2  X  Rchip  =  2.046  MHz. 

The  prompt  eode  Cp{t)  is  generated  using  the  same  spreading  eode  sequenee  used  by 
the  DSSS  transmitter  but  with  its  timing  eharaeteristies  varied  by  the  NCO  in  aeeordanee 
with  the  DLL  error  voltage.  The  early  eode  Ce{t)  is  advaneed  Tchip/2  relative  to  the  prompt 
eode  while  the  late  eode  ci{t)  is  delayed  Tchip/2  relative  to  the  prompt  eode. 

The  IF  bandpass  filters  serves  to  remove  all  but  the  fundamental  harmonie  after  mul- 
tiplieation  by  either  Cc{t)  or  ci{t).  The  IF  filters  were  implemented  as  8*^-order  Butterworth 
filters  having  a  —3  dB  bandwidth  of  Wbpf  =  5  kHz.  This  bandwidth  is  somewhat  wider 
than  neeessary  to  traek  the  eode  of  interest,  but  a  wider  bandpass  was  ehosen  for  eonsis- 
teney  with  real-world  eonditions  where  Doppler  shift  is  not  preeisely  known  or  not  traeked. 

The  squaring  operation  following  the  IF  bandpass  filtering  and  subsequent  low  pass 
filtering  eomprise  an  envelope  deteetor.  The  low  pass  filters  are  designed  to  remove  dou¬ 
ble  frequeney  terms  resulting  from  the  squaring  operation  and  were  implemented  here  4*^- 
order  Chebyehev  filters  having  a  —3  dB  bandwidth  of  Wlpf  =  2.5  MHz.  The  early  and  late 
gate  low  pass  filter  outputs  are  summed  and  filtered  to  provide  the  eontrol  signal  (diserim- 
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Figure  4.7:  Delay-Lock  Loop  (DLL)  architecture  used  for  SMSE-DSSS 
coexistence  demonstration  under  imperfect  code  tracking  conditions. 


inator  voltage)  to  the  Numerically  Controlled  Oscillator  (NCO).  For  simulation  purposes, 
the  NCO  maintains  its  own  time  reference  to  generate  the  PN  code.  When  the  discrimina¬ 
tor  voltage  is  zero,  the  NCO  sample  rate  matches  the  simulated  sample  rate  of  fs  and  the 
estimated  prompt  code  Cp{t)  perfectly  matches  the  received  code.  When  the  signed  discrim¬ 
inator  value  deviates  positive/negative  from  zero,  the  NCO  sample  rate  increases/decreases 
accordingly  and  the  estimated  prompt  code  Cp{t)  is  altered.  Ideally,  the  discriminator  value 
varies  linearly  between  its  extremes.  When  viewed  in  real-time,  an  accordion-like  effect  is 
observed  in  Cp{t)  on  a  chip-by-chip  basis,  i.e.,  there  is  a  clear  compression  and  dilation  of 
the  chip  intervals  throughout  the  code  period. 

Perfect  DLL  tracking  results  in  a  prompt  code  estimate  that  is  an  exact  replica  (code 
phase  and  chip  duration)  of  the  received  DSSS  spreading  waveform.  The  presence  of  a 
coexisting  SMSE  signal  causes  the  estimated  prompt  code  to  vary  in  both  phase  and  chip 
duration  as  the  DLL  compensates.  This  degradation  is  characterized  and  quantified  using 
the  cross-correlation  metric  Coeg  introduced  in  Section  3.2.3. 

The  correlation  metric  introduced  in  Section  3.2.3  is  generated  as  follows  to  charac¬ 
terize  DLL  code  tracking  performance.  The  DSSS  waveform  is  received  by  the  DLL  with 
an  initial  code  phase  difference  of  20%  between  the  prompt  code  and  the  true  spreading 
code.  The  DLL  processes  this  signal  for  two  code  cycles  without  adaptation  to  allow  the 
DLL  to  stabilize  to  an  accurate  initial  phase  estimate.  Adaption  begins  after  two  code  pe¬ 
riods.  After  five  code  periods  the  coexistent  SMSE  signal  is  applied  to  the  loop  along  with 
the  DSSS  signal  for  five  additional  code  periods.  The  prompt  code  and  true  spreading  code 
from  the  five  additional  code  periods  are  correlated  per  (3.2)  using  an  integration  interval 
of  one  code  period.  The  mean  value  from  the  correlation  process  comprises  one  sample  of 
the  correlation  metric,  Ct-p-  Eor  clarity,  the  data  is  presented  is  terms  of  a  normalized  and 
shifted  version  of  Cp-p,  termed  Cpeg^  which  represents  the  amount  of  tracking  degradation 
caused  by  the  coexistent  SMSE  signal. 

Representative  histograms  of  the  Cpeg  metric  are  shown  in  Eig.  3.1  and  Eig.  3.2  for 
moderately  and  severely  degraded  DEE  code  tracking  performance,  respectively.  Recall 
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that  values  of  Coeg  approaching  zero  indicate  better  DLL  tracking  performance  while  val¬ 
ues  near  one  indicate  poorer,  more  degraded  tracking  performance. 

In  addition,  observations  of  Coeg  behavior  in  many  such  figures  indicates  that  the 
SMSE  waveform  impact  on  DLL  performance  is  highly  dependent  on  ISR.  Thus,  care  was 
taken  to  choose  appropriate  ISR  values  for  subsequent  waveform  design  demonstrations. 
For  example,  if  the  chosen  ISR  value  is  too  high  all  of  the  potential  SMSE  waveform 
designs  can  cause  the  DEE  to  break  lock.  Eikewise,  if  the  chosen  ISR  value  is  too  low, 
many  of  the  potential  waveform  designs  will  have  minimal  impact.  Observations  of  C^eg 
for  15  dB  <  ISR  <  35  dB  showed  that  the  metric  behaved  best  for  optimization  purposes 
at  ISR  =  20  dB.  At  that  value,  all  of  the  SMSE  parameter  combinations  degraded  DEE 
code  tracking  performance,  but  none  consistently  caused  the  DEE  to  break  lock. 

4.2.2  Hybrid  Optimization  Step  1:  GA  Process.  The  hybrid  optimization  tech¬ 
nique  first  uses  the  GA  process  in  Section  3.2.4  to  find  a  “coarse”  optimized  SMSE  solution. 
For  demonstration  purposes,  1000  generations  were  used  with  optimization  variables  being 
the  number  of  IFFT  points  Nf  G  {16,  32,  64, 128}  and  the  subcarrier  frequency  spacing 
A/  e  [1,  33]  n  X.  The  feasibility  region  was  defined  such  that  NfAf  <  1/Tchip- 

The  GA  process  was  used  to  both  minimize  and  maximize  G^eg  under  imperfect  DEE 
tracking  conditions.  Table  4.6  shows  the  initial  random  and  final  optimized  populations 
using  the  GA  process  to  minimize  Coeg  ■  Box  and  whisker  plots  of  Coeg  for  the  initial  and 
final  optimized  populations  are  shown  in  Fig.  4.8  and  Fig.  4.9,  respectively. 

Table  4.7  shows  the  initial  random  and  final  optimized  populations  using  the  GA 
process  to  maximize  Coeg-  Box  and  whisker  plots  of  Coeg  for  the  initial  and  final  optimized 
populations  are  shown  in  Fig.  4.10  and  Fig.  4.11,  respectively. 

The  parameter  combinations  in  Table  4.6  and  Table  4.7  were  compared  via  the  ESD 
test  described  in  Section  2. 2.2. 8.  The  SMSE  parameters,  (Nf  =  16,  A/  =  6  kHz),  pro¬ 
vided  Coeg  means  which  were  statistically  lower  than  all  other  parameter  combinations  ex¬ 
cept  for  (Nf  =  16,  A/  =  7  kHz).  These  two  SMSE  parameter  combinations  were  statisti¬ 
cally  indistinguishable.  The  combination  (Nf  =  16,  A/  =  6  kHz)  was  chosen  as  the  start- 


71 


Table  4.6:  Initial  random  and  final  optimized  SMSE  parametrie  value 
populations  using  the  GA  proeess  to  MINIMIZE  eorrelation  degradation 
Coeg  -  Best-Case  Traeking  and  SMSE-DSSS  Coexistenee. 


Nf 

A/  {kHz) 

Eigs.  4.8  &  4.9 

Initial 

Optimized 

Initial 

Optimized 

Combination 

32 

16 

23 

10 

1 

64 

16 

8 

38 

2 

16 

32 

53 

6 

3 

16 

16 

17 

6 

4 

32 

16 

14 

6 

5 

16 

32 

24 

6 

6 

64 

16 

11 

6 

7 

16 

16 

7 

6 

8 

16 

16 

34 

6 

9 

16 

16 

63 

6 

10 

Table  4.7:  Initial  random  and  final  optimized  SMSE  parametrie  value 
populations  using  the  GA  proeess  to  MAXIMIZE  eorrelation  degradation 
Coeg  -  Worst-Case  Traeking  and  SMSE-DSSS  Coexistenee. 


Nf 

A/  {kHz) 

Eigs.  4.10  &  4.11 

Initial 

Optimized 

Initial 

Optimized 

Combination 

16 

32 

48 

25 

1 

16 

32 

20 

24 

2 

16 

32 

24 

28 

3 

64 

32 

13 

24 

4 

32 

32 

29 

24 

5 

32 

32 

27 

24 

6 

32 

32 

31 

24 

7 

16 

32 

33 

24 

8 

32 

32 

22 

24 

9 

16 

32 

29 

24 

10 
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Figure  4.8:  Box  and  whisker  plot  of  Coeg  results  for  initial  random  pop 
ulation  using  the  GA  process  to  MINIMIZE  Coeg- 
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Figure  4.9:  Box  and  whisker  plot  of  Coeg  results  fox  final  optimized  pop 
ulation  using  the  GA  process  to  MINIMIZE  Coeg- 
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Table  4.8:  Hybrid  Optimization  Process:  GA  Optimized  SMSE  param¬ 
eters  Nf  and  A/  maximizing  and  minimizing  the  correlation  degradation 
metric,  Coeg- 


Nf 

A/  {kHz) 

Wsmse  {MHz) 

Minimized  Coeg 

16 

6 

0.19 

Maximized  G^eg 

32 

24 

1.54 

ing  point  for  the  RSM  minimization  process  instead  of  (A/  =  16,  A/  =  7  kHz),  because  it 
was  the  result  of  GA  convergence.  However,  the  parameter  values  (A/  =  16,  A/  =  7  kHz) 
were  also  explored  with  the  RSM  process,  because  of  its  proximity  to  (Aj  =  16,  A/  = 
6  kHz). 

The  SMSE  parameters  (A/  =  32,  A/  =  24  kHz)  provided  Coeg  means  that  were 
not  statistically  lower  than  any  other  parameter  combinations.  Eleven  SMSE  parame¬ 
ter  combinations  from  Table  4.6  and  Table  4.7  provided  means  which  were  statis¬ 
tically  indistinguishable  from  that  of  (A/  =  32,  A/  =  24  kHz).  The  other  parameter 
combinations  were  similar  to  {Nf  =  32,  A/  =  24  kHz)  with  respect  to  their  relatively 
high  bandwidths  (Wsmse  =  ^  ^/)-  The  bandwidth  range  for  these  signals  was 

512kHz  <  Wsmse  <  992  kHz.  Eor  comparison  purposes,  the  bandwidth  of  the  mini¬ 
mization  solution,  {Nf  =  16,  A/  =  6  kHz)  was  Wsmse  =  96  kHz.  The  combination 
{Nf  =  32,  A/  =  24  kHz)  was  chosen  as  the  starting  point  for  the  RSM  maximization 
process  instead  of  the  other,  similarly  performing  SMSE  parameter  combinations,  because 
it  was  the  result  of  GA  convergence. 

The  resultant  GA-optimized  SMSE  parameters  are  shown  in  Table  4.8,  and  their 
corresponding  histograms  for  Coeg  are  shown  in  Eig.  4.12  and  Eig.  4.13. 

4.2.3  Hybrid  Optimization  Step  2:  RSM  Process.  Given  the  “coarse”  optimiza¬ 
tion  solution  from  the  GA  process  in  Section  4.2.2,  the  RSM  process  in  Section  3.2.5  is 
next  applied  to  perform  final  minimization  and  maximization  of  correlation  degradation 
Coeg-  This  hybrid  GA-RSM  approach  addressed  two  issues  that  emerged  in  Section  4.1.3. 
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Figure  4.10:  Box  and  whisker  plot  of  Coeg  results  for  initial  random 
population  using  the  GA  proeess  to  MAXIMIZE  Coeg- 
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Figure  4.11:  Box  and  whisker  plot  of  Coeg  results  for  final  optimized 
population  using  the  GA  proeess  to  MAXIMIZE  Cneg- 
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Correlation  Degragation,  Coeg 


Figure  4.12:  Histogram  showing  behavior  of  Coeg  for  severely  degraded 
DLL  code  tracking  with  ISR  =  25  dB  (at  the  output  of  the  front-end  RF 
filter)  in  coexistent  SMSE-DSSS  scenario  with  GA  MAXIMIZED  SMSE 
Parameters  {Nf,Af). 


Correlation  Degragation,  Coeg 


Eigure  4.13:  Histogram  showing  behavior  of  C^eg  for  moderately  de¬ 
graded  DEE  code  tracking  with  ISR  =  25  dB  (at  the  output  of  the  front-end 
RE  filter)  in  coexistent  SMSE-DSSS  scenario  with  GA  MINIMIZED  SMSE 
Parameters  {Nf,Af). 
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Table  4.9:  Hybrid  Optimization  Process:  SMSE  parameter-to-coded  vari¬ 
able  (A/-to-a:)  mapping  for  initial  RSM  MINIMIZATION  experiment  under 
imperfect  code  tracking  conditions. 


X 

A/  {kHz) 

-2 

6.00 

-1 

6.23 

0 

6.47 

1 

6.71 

2 

6.94 

The  first  issue  had  to  do  with  answering  the  question,  “What  is  a  good  starting  point  for 
the  RSM  process?”  Given  GA  is  the  first  step  in  the  hybrid  process,  a  systematic  approach 
is  now  used  to  choose  the  RSM  starting  point.  The  second  issue  is  simply  related  to  the 
number  of  dimensions  used  in  the  optimization  process.  Given  the  GA  solution  effectively 
identifies  the  appropriate  Nj  value,  the  RSM  optimization  process  is  reduced  to  a  one  di¬ 
mensional  search  for  the  appropriate  A/  value.  This  is  a  desirable  outcome  given  that  RSM 
attempts  to  optimize  on  a  continuous  surface,  while  N f  is  inherently  discrete.  Thus,  resul¬ 
tant  RSM  design  matrix  has  one  factor  with  five  levels.  The  number  of  SMSE  subcarriers 
was  set  to  Aj  =  16  as  determined  by  the  initial  GA  process. 

4.2.3. 1  Minimizing  DLL  Tracking  Degradation.  Einal  RSM  optimization 
is  first  considered  for  the  case  where  Coeg  is  to  be  minimized  for  the  coexistent  SMSE- 
DSSS  scenario.  The  values  used  for  A/  are  shown  in  Table  4.9  along  with  their  coded 
counterparts  in  variable  x. 

The  RSM  design  matrix  was  run  with  1000  repetitions  and  the  data  transformed  using 
the  Box-Cox  transformation  with  a  value  of  A  =  54.78.  This  resulted  in  each  condition  in 
the  trial  being  properly  classified  as  normal  according  to  the  Shapiro-Wilks  test  described 
in  Section  2. 2. 2. 5. 

The  data  was  then  fit  to  linear,  pure  quadratic,  and  quadratic  models  as  described  in 
Section  2. 2.2. 2.  The  resultant  p-value  test  for  all  three  models  yielded  p  <  10“^  indicating 
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Table  4.10:  Hybrid  Optimization  Process:  SMSE  parameter-to-coded 

variable  (Af-to-x)  mapping  for  initial  RSM  MAXIMIZATION  experiment 
under  imperfect  code  tracking  conditions. 


X 

A/  {kHz) 

-2 

23.03 

-1 

23.52 

0 

24.00 

1 

24.48 

2 

24.97 

a  good  fit.  The  pure  quadratic  model  provided  the  best  fit  based  on  its  F  statistic  and 
p- value.  However,  the  solutions  for  all  three  models  were  calculated  since  all  were  deemed 
to  be  statistically  significant. 

Based  on  these  results,  the  final  solution  for  a  coexistent  SMSE  signal  that  is  least  dis¬ 
ruptive  to  DSSS  system  performance,  as  indicated  by  minimal  DSSS  code  tracking  degra¬ 
dation,  is  Nf  =  16  and  A/  =  6.47  kHz.  This  was  accomplished  by  finding  a  coarse  GA 
solution  followed  by  fine-tuning  with  the  RSM  process. 

4. 2. 3. 2  Maximizing  DLL  Tracking  Degradation.  Einal  RSM  optimization 
is  first  considered  for  the  case  where  Coeg  is  to  be  maximized  for  the  coexistent  SMSE- 
DSSS  scenario.  The  values  used  for  A/  are  shown  in  Table  4.10  along  with  their  coded 
counterparts  in  variable  x. 

The  RSM  design  matrix  was  run  with  1000  repetitions  and  the  data  transformed  using 
the  Box-Cox  transformation  with  a  value  of  A  =  10.88.  This  resulted  in  each  condition  in 
the  trial  as  being  properly  classified  as  normal  according  to  the  Shapiro-Wilks  test  described 
in  Section  2. 2. 2. 5. 

The  data  was  then  fit  to  linear,  pure  quadratic,  and  quadratic  models  as  described 
in  Section  2. 2. 2. 2.  The  resultant  p-value  test  for  the  quadratic  model  yielded  p  <  10“^ 
indicating  a  good  fit. 
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Table  4.11:  Hybrid  Optimization  Process:  RSM  Optimized  SMSE  pa¬ 
rameters  Nf  and  A/  maximizing  and  minimizing  the  correlation  degrada¬ 
tion  metric,  Coeg- 


Nf 

Af  (kHz) 

Wsmse  {MHz) 

Minimized  Coeg 

16 

6.47 

0.2 

Maximized  Cneg 

32 

25.02 

1.6 

Based  on  these  results,  the  final  solution  for  a  coexistent  SMSE  signal  that  is  the 
most  disruptive  to  DSSS  system  performance,  as  indicated  by  maximal  DSSS  code  track¬ 
ing  degradation,  is  Nf  =  32  and  A/  =  25.02  kHz,  as  shown  in  Table  4.11.  This  was 
accomplished  by  finding  a  coarse  GA  solution  followed  by  fine-tuning  with  the  RSM  pro¬ 
cess.  To  visually  compare  the  final  solutions  of  the  RSM  process,  histograms  of  the  final 
minimized  and  maximized  Cneg  values  are  provided  in  Eig.  4.14  and  Eig.  4.15,  respectively. 

Given  the  resultant  Wsmse  for  each  optimized  case  in  Table  4.11,  all  SMSE  power 
is  contained  within  Wrf  =  2.046  MHz  and  thus  the  RE  filtering  effects  contribute  min¬ 
imally  to  the  performance  difference.  The  performance  difference  is  most  attributable  to 
the  relationship  between  SMSE  subcarrier  spacing  and  the  spectral  structure  of  the  DSSS 
spreading  signal,  i.e.,  the  degree  of  SMSE-DSSS  spectral  coincidence.  Eor  parameters 
used  in  these  imperfect  tracking  scenarios,  the  spectral  characteristics  of  the  DSSS  spread¬ 
ing  signal  and  received  SMSE  signals  can  be  summarized  as  follows: 


1.  The  spectral  lines  for  the  Nc  =  1023-chip  DSSS  spreading  code  are  spaced  A  f chip  = 
Rchip/Nc  =  1  KHz  apart  [25].  There  are  Wrf/ Af chip  =  2,  046  total  lines  in  Wrf  = 
2  X  Rchip  =  2.046  MHz  that  are  power  weighted  according  to  [25] 


(/  -  fc) 


1 

ATI’ 


RffA  sinc^  [{f  -  fc)  Tchip] , 


f  =  fc 

Elsewhere 


where  sinc(/)  =  sin  {f)/ f  ■  The  central  line  is  located  at  fc  and  remaining  lines 
uniformly  spaced  on  either  side  of  fc  at  intervals  of  Afchip  =  1  KHz. 
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Correlation  Degragation,  Coeg 


Figure  4.14:  Histogram  showing  behavior  of  Coeg  for  the  RSM  solution 
that  MINIMIZES  Coeg- 
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Correlation  Degragation,  Coeg 

Figure  4.15:  Histogram  showing  behavior  of  Coeg  for  the  RSM  solution 
that  MAXIMIZES  Cueg- 
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2.  The  coexisting  SMSE  signal  spectrum  is  centered  at  fc  with  an  even  number  of  BPSK 
data  modulated  subcarriers  (Nf/2)  uniformly  spaced  on  either  side  of  fc-  The  two 
centrally  located  SMSE  subcarriers  are  at  A/  on  either  side  of  fc  and  all  remaining 
subcarriers  spaced  at  uniform  intervals  of  A/  from  these.  By  design,  the  total  SMSE 
power  /  is  equally  distributed  across  all  subcarriers. 

The  degree  of  spectral  coincidence  between  the  DSSS  spreading  code,  as  used  in  the 
DSSS  receiver  to  “despread”  the  received  signals,  and  received  SMSE  signals  dictates  Pb 
performance.  This  is  best  characterized  by  considering  the  relationship  between  the  DSSS 
spectral  line  spacing  Afchip  and  the  SMSE  subcarrier  spacing  A/,  or  more  directly,  by 
considering  the  resultant  power  spectrum  at  the  despreading  mixer  output  (DSSS  detector 
input)  in  the  DSSS  receiver.  The  “despread”  spectral  response  is  the  convolution  of  the 
uniformly  spaced  DSSS  spectral  lines  at  Afchip  =  1  KHz  with  the  uniformly  spaced 
SMSE  subcarriers  spaced  at  A/. 

Eor  the  least  disruptive  solution  in  Table  4.11  the  optimized  subcarrier  spacing  of 
A/  =  6.47  KHz  is  a  non-integer  multiple  of  Afchip  =  1  KHz.  Therefore,  for  a  given 
shift-multiply-integrate  operation  of  the  spectral  convolution  process  there  is  a  most  one 
DSSS  spectral  line  that  is  aligned  with  an  SMSE  subcarrier  and  the  peak  response  of  all 
other  SMSE  subcarriers  falls  approximately  midway  between  other  spectral  lines-minimal 
SMSE-DSSS  spectral  coincidence.  As  a  result,  the  amount  of  power  in  the  “despread” 
spectral  response  that  falls  within  the  DSSS  detector  bandwidth  (W oet  =  2  x  Rsym  = 
2  KHz)  is  relatively  low  and  the  interfering  SMSE  signal  has  minimal  impact  on  DSSS  Pb 
performance. 

On  the  other  hand,  for  the  most  disruptive  solution  in  Table  4.1 1,  the  optimized  sub¬ 
carrier  spacing  of  A/  =  25.02  KHz  is  approximately  an  integer  multiple  of  Afchip  = 
1  KHz.  Therefore,  for  a  given  shift-multiply-integrate  operation  of  the  spectral  convolu¬ 
tion  process  the  peak  responses  of  all  SMSE  subcarriers  align  near  perfectly  with  specific 
DSSS  spectral  lines-maximum  SMSE-DSSS  spectral  coincidence.  As  a  result,  the  amount 
of  power  in  the  “despread”  spectral  response  that  falls  within  the  DSSS  detector  bandwidth 
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(Woet  =  2  X  Rsym  =  2  K H z)  IS  maximum  and  the  interfering  SMSE  signal  has  maxi¬ 
mum  impaet  on  DSSS  Pb  performanee.  Alternately  stated,  this  SMSE  solution  in  this  ease 
is  more  speetrally  “matehed”  to  the  DSSS  signal  which  causes  greater  degradation. 

4.2.4  Optimal  SMSE  Parameter  Demonstration.  Two  additional  metrics  were 
calculated  to  demonstrate  waveform  optimization  effects  using  the  final  RSM  optimized 
solutions  from  Section  4.2.  Eirst,  DSSS  end-to-end  bit  error  rate  Pb  was  calculated  to 
demonstrate  the  communication  system  applicability.  Second,  relative  timing  jitter,  Jcode, 
between  the  DEE- generated  prompt  code  Cp{t)  and  the  incident  DSSS  code  Ci{t)  was  calcu¬ 
lated  to  demonstrate  applicability  to  DSSS-based  navigation  and  precision  timing  system. 

Communication  system  Pb  results  are  shown  in  Table  4.12.  These  results  were  gen¬ 
erated  using  the  physics-based  analytic  DSSS  receiver  model  with  Monte  Carlo  simulation 
as  shown  in  Section  4.1.  The  ideal  results  are  provided  for  comparison  and  represent  per¬ 
fect  code  tracking  conditions,  i.e.,  the  best  that  can  be  expected.  The  minimization  and 
maximization  results  were  obtained  using  the  SMSE  parameter  solutions  from  the  hybrid 
optimization  process.  Eor  the  degraded  cases,  the  SMSE  waveform  effects  on  the  DSSS 
receiver  were  isolated  to  DEE  tracking  performance  only;  no  SMSE  signal  was  present  in 
the  DSSS  demodulator. 

Timing  jitter  Jcode  results  are  shown  in  Table  4.13.  Eor  this  analysis,  Jcode  was  defined 
as  the  the  time  difference  between  the  coded  pulse  transition  points  in  Ci{t)  and  Cp{t),  with 
negative  values  for  Jcode  indicating  that  Cp{t)  transitioned  before  Cp{t)  and  positive  values 
indicating  that  Cp{t)  transitioned  after  Cp{t).  Table  4.13  shows  the  standard  deviation  of  the 
Jcode  metric. 
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Table  4.12:  Communication  system  bit  error  rate  (Pfe)  for  ideal  code  track¬ 
ing  and  degraded  code  tracking  using  final  SMSE  parameters  from  hybrid 
GA-RSM  optimization  process  in  Section  4.2. 


Coeg  Optimization 

Nf 

AfikHz) 

Pb 

Ideal 

N/A 

N/A 

1.0  X  10-3 

Minimized 

16 

6.47 

5.6  X  10-2 

Maximized 

32 

25.02 

2.2  X  10-^ 

Table  4.13:  Timing  jitter  standard  deviation  Jcode  for  degraded  code  track¬ 
ing  using  final  SMSE  parameters  from  hybrid  GA-RSM  optimization  pro¬ 
cess  in  Section  4.2. 


Coeg  Optimization 

Nj 

AfikHz) 

Jcode  (56c) 

Minimized 

Maximized 

16 

32 

6A1 

25.02 

1.83  X  10-^ 
1.94  X  10-^ 
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V.  Conclusion 


5.1  Research  Summary 

The  Spectrally  Modulated,  Spectrally  Encoded  (SMSE)  framework  provides  an  ef¬ 
fective  means  for  implementing  Orthogonal  Erequency  Division  Multiplexing  (OEDM)  sig¬ 
nals  and  the  ability  to  efficiently  generate  them  has  only  recently  become  practical  from  a 
hardware  perspective.  As  such,  OEDM  technology  has  been  identified  as  a  bedrock  tech¬ 
nology  for  fourth  generation  (4G)  communications  based  on  Cognitive  Radio  (CR)  and 
Software  Defined  Radio  (SDR)  techniques  [17].  The  inherent  flexibility  of  OEDM  has 
generated  significant  research  interest  [1 1,20,34-37,40]  and  the  expanding  pool  of  OEDM 
variants  drove  the  need  for  a  unified  framework.  As  utilized  under  this  research,  the  resul¬ 
tant  SMSE  framework  effectively  embodies  the  class  of  OEDM-based  signals  [27-30]. 

As  4G  SMSE  communications  emerge  they  must  coexist  with  other  systems  while 
competing  to  use  available  communication  resources.  Thus,  coexistent  interference  be¬ 
comes  a  concern,  especially  when  these  systems  operate  without  orthogonality  in  fre¬ 
quency,  time,  space,  polarization,  and/or  coding.  Given  a  lack  of  orthogonality  and  limited 
available  bandwidth,  these  signals  must  be  designed  to  spectrally  coincide  while  induc¬ 
ing  “manageable”  levels  of  mutual  interference.  This  becomes  particularly  challenging 
when  fundamentally  different  waveform  modulations  and  the  waveform  design  procedure 
often  resorts  to  trial  and  error  design  methods  [25,  31].  The  goal  of  this  research  was  to 
demonstrate  a  more  structured,  optimal  means  for  SMSE  waveform  design  using  tech¬ 
niques  commonly  employed  in  the  operations  research  field.  A  survey  of  general  optimiza¬ 
tion  techniques  revealed  that  two  methods  were  particularly  applicable  to  the  coexistent 
SMSE  waveform  design  scenario,  including  Genetic  Algorithm  (GA)  and  Response  Sur¬ 
face  Methodology  (RSM)  optimization  techniques. 

Each  of  these  optimization  techniques  are  used  to  demonstrate  SMSE  waveform  de¬ 
sign  in  a  coexistent  scenario  containing  an  SMSE  signal  and  a  Direct  Sequence  Spread 
Spectrum  (DSSS)  system  operating  over  an  Additive  White  Gaussian  Noise  (AWGN)  chan¬ 
nel.  The  DSSS  system  employs  a  non-coherent  Delay-Eock  Eoop  (DEE)  for  code  track¬ 
ing.  The  specific  DEE  implementation  used  here,  along  with  defined  correlation  metrics 
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that  characterize  eode  traeking  eonditions,  is  suffieiently  general  sueh  that  the  optimization 
demonstrations  herein  are  broadly  applieable  to  other  non-eommunieation  applieations  em¬ 
ploying  DLL  traeking,  e.g.,  preeision  navigation,  timing,  geoloeation,  ete. 

The  eoexistent  SMSE-DSSS  waveform  design  proeess  is  demonstrated  herein  under 
both  perfect  and  imperfect  DLL  eode  traeking  eonditions.  Under  both  eonditions  the  num¬ 
ber  of  subcarriers  N f  and  subcarrier  spaeing  A/  are  the  SMSE  design  variables  of  interest. 
Under  perfect  DEE  eode  traeking  eonditions,  the  GA  and  RSM  optimization  proeesses 
are  eonsidered  independently  with  the  objeetive  funetion  to  be  optimized  (minimized  and 
maximized)  being  DSSS  bit  error  rate  P^.  A  hybrid  GA-RSM  optimization  proeess  is  used 
under  more  realistie  imperfect  DEE  eode  traeking  eonditions.  In  this  ease,  Pb  optimization 
is  aeeomplished  through  a  eorrelation  degradation  metrie  with  the  GA  proeess  being  first 
applied  to  generate  a  “eoarse”  solution  followed  by  RSM  proeessing  to  provide  the  final 
optimized  solution. 

Eor  all  perfeet  and  imperfeet  DEE  eode  traeking  seenarios  eonsidered,  the  optimized 
DSSS  Pb  minimization  results  yielded  SMSE  waveform  designs  and  Pb  performanee  that 
was  eonsistent  with  seenarios  having  no  eoexistent  SMSE  signal  present  (best-ease  eo¬ 
existent  performanee).  Eor  the  optimized  DSSS  Pb  maximization  solutions,  worst-ease 
SMSE-DSSS  eoexistenee  was  aehieved  for  SMSE  waveform  designs  that  were  speetrally 
“matehed”  to  the  DSSS  signal,  i.e.,  greatest  Pb  degradation  was  experieneed  when  the  re¬ 
sultant  SMSE  subearrier  spaeing  A/  was  an  integer  multiple  of  the  speetral  line  spaeing 
^fchip  of  the  DSSS  spreading  eode. 

The  researeh  objeetive  has  been  aehieved  in  the  sense  that  4G  eommunieations  de¬ 
sign  engineers  now  have  one  additional  tool  at  their  disposal.  This  work  has  suoeessfully 
expanded  the  praetieal  utility  of  a  previously  developed  tool,  the  original  SMSE  frame¬ 
work  [26,28,30],  by  demonstrating  a  more  effieient,  struetured  means  for  eoexistent  wave¬ 
form  design  that  replaees  previous  trial  and  error  methods.  As  sueh,  the  eommunieations 
eommunity  is  one  step  eloser  to  aetually  hitting  the  bedroek  of  OEDM-based  signaling 
using  the  SMSE  framework.  The  signifieanee  of  this  has  been  aeknowledged  through  ae- 
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ceptance  of  the  independent  GA  [3]  and  RSM  [4]  optimization  results  under  perfect  DLL 
eode  tracking  conditions.  Applicability  of  the  hybrid  GA-RSM  processing  technique  under 
imperfect  DLL  code  tracking  conditions  has  been  acknowledged  as  well  [2] . 

5.2  Recommendations  for  Future  Research 

Given  that  applicability  of  optimization  techniques  to  SMSE  waveform  design  has 
been  demonstrated  and  well-received  by  the  technical  community,  there  are  many  addi¬ 
tional  research  topics  that  could  be  investigated.  Some  of  the  most  evident  future  research 
avenues  include: 

•  The  research  assumptions  in  Chapter  I  could  be  removed  and/or  relaxed  and  the 
research  process  repeated.  In  some  cases,  this  would  represent  a  somewhat  triv¬ 
ial  exercise  and  produce  results  that  might  be  expected.  For  example,  the  coexis¬ 
tent  SMSE-DSSS  demonstration  results  in  Chapter  V  could  easily  be  expanded  by 
considering  alternative  OFDM-based  communication  variants  (CI-OFDM,  COFDM, 
MC-CDMA,  etc.).  In  other  cases,  relaxing  the  research  assumptions  could  prove  to 
be  far  more  challenging  and  could  produce  results  which  are  far  more  significant. 

•  The  demonstrations  here  focused  on  SMSE  waveform  design  through  appropriate 
selection  of  two  variables,  Nj  and  A/.  For  all  cases,  the  total  SMSE  power  was 
equally  distributed  across  all  selected  subcarriers.  Research  could  be  conducted  that 
maintains  the  same  coexistent  SMSE-DSSS  scenarios  and  GA/RSM  optimization 
objective  functions  but  with  additional  consideration  given  to  incorporating  the  se¬ 
lection  of  SMSE  weight  vector  Wk.  By  treating  each  element  of  Wk  as  a  model 
parameter,  the  ANOVA  process  could  provide  insight  into  the  impact  that  individual 
subcarriers  are  having  on  the  objective  function.  Assuming  the  optimization  goal  is  to 
achieve  best-case  coexistence,  the  practical  implication  is  that  insignificant  subcarri¬ 
ers  (those  inducing  minimal  interference)  could  transmit  at  appreciable  power  levels 
while  power  levels  in  significant  subcarriers  (those  inducing  maximum  interference) 
could  be  reduced  or  set  to  zero  (shut  off). 
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•  The  approach  to  shutting  off  significant  subcarriers  is  perhaps  best  characterized  as 
being  a  hard  decision  process.  Recent  efforts  have  introduced  overlay,  underlay 
and  hybrid  overlay /underlay  waveform  design  using  the  original  SMSE  framework 
with  subcarrier  selection  and  power  distribution  based  on  soft  decision  criteria  [5, 6]. 
Given  the  inherent  structure  of  the  original  SMSE  framework  is  maintained  in  these 
efforts,  the  emergence  of  these  new  SD-SMSE  techniques  warrants  future  investi¬ 
gation  given  that  the  optimization  methods  demonstrated  in  this  work  should  to  be 
directly  applicable. 

•  Demonstrations  in  this  work  are  based  on  internally  observable  knowledge  of  how 
the  coexistent  DSSS  system,  and  in  particular  its  DEE  code  tracking,  is  respond¬ 
ing  to  various  SMSE  parameter  changes.  This  internal  knowledge  will  generally 
be  unavailable  and  thus  alternatives  need  to  be  investigated.  Two  potential  alterna¬ 
tives  could  be  considered.  Eirst,  cooperative  exchange  of  system  performance  status 
among  coexistent  system(s)  through  communication  back  channels.  While  this  ap¬ 
proach  generally  requires  additional  communication  resources,  which  may  result  in 
fewer  resources  being  available  for  the  primary  function,  similar  techniques  have 
been  employed  in  communication  networks.  Second,  the  objective  function  behav¬ 
ior  could  be  determined  passively  without  any  aid  from  coexistent  systems.  While 
this  approach  is  technically  more  challenging,  it  is  consistent  with  envisioned  4G 
communication  goals  for  CR-based  SDR  communications  that  will  use  externally 
observable  knowledge  to  adapt  and  optimize  performance. 

•  The  specific  non-coherent  DEE  implementation  used  here,  along  with  defined  cor¬ 
relation  metrics  for  characterizing  code  tracking  conditions,  are  sufficiently  gen¬ 
eral  such  that  the  optimization  processes  considered  are  broadly  applicable  to  other 
non-communication  applications  employing  DEE  tracking,  e.g.,  precision  naviga¬ 
tion,  timing,  geolocation,  etc.  Therefore,  a  similar  SMSE  coexistence  analysis  and 
demonstration  could  be  conducted  within  each  of  these  application  areas. 
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